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DIAGNOSIS  OF  THE  COMMA  CLOUD 
OF  10  APRIL  1979 

CHAPTER  I 

INTRODUCTION 

The  cloud  patterns  associated  with  synoptic  scale,  midlatitude 
cyclonic  disturbances  usually  take  on  a  spiral  shape  equatorward  of  the 
upper  cyclone  with  a  broader  region  of  cloudiness  spreading  poleward 
across  the  path  of  the  upper  trough.  Often  the  pattern  is  isolated  from 
other  cloudy  regions,  taking  the  shape  of  a  vast  comma  ~2000km  long;  and 
by  common  usage,  any  vaguely  similar  cloud  pattern  that  accompanies 
short  waves  in  the  westerlies  has  come  to  be  known  as  a  "comma  cloud." 

Unfortunately,  the  term  "comma  cloud"  is  used  also  to  refer  to 
meso-oc cloud  vortices  -600km  long,  which  occur  with  vorticity  maxima  in 
the  relatively  cold  air  around  upper  level  extratropical  cyclones. 

Boucher  and  Newcomb  (1962)  proposed  a  model  of  the  lifecycle 
of  the  synoptic  scale  vortex  patterns  related  to  the  evolving  surface 
frontal  wave  with  which  they  are  associated.  Efforts  to  infer  the  ver¬ 
tical  motion  field  in  the  vicinity  of  comma  clouds  directly  from 
satellite  observed  cloudiness  failed  (Timchalk  and  Hubert  (1961),  Leese 
(1962),  Hansen  and  Thompson  (1964),  Nagle,  et  al .  (1966),  Barr,  et  al . 


(1966)),  but  the  work  demonstrated  that: 

(1)  Instantaneous  vertical  velocity  fields  do  not  agree  with 
synoptic-scale  cloud  patterns,  except  possibly  In  the  early  stages  of 
storm  development.  Maxima  of  instantaneous  upward  and  downward  motion 
tend  to  be  upstream  of  the  centers  of  cloudy  and  cloud-free  areas, 
respectively,  especially  in  rapidly  moving  systems.  This  reinforces  the 
view  that  clouds  result  from  the  time-integrated  ascent  of  parcels  which 
become  saturated  somewhere  downstream  of  the  maximum  of  instantaneous 
ascent. 

(2)  Other  factors,  such  as  moisture  stratification  and  the 
history  of  vertical  motion,  which  air  parcels  have  experienced,  control 
the  development  and  dissipation  of  clouds. 

(3)  The  shape  of  major  cloud  patterns  is  due  mainly  to  the 
horizontal  advection  of  pre-existing  cloud  matter. 

This  study  will  attempt  to  verify,  as  did  Barr,  et  al .  (1966), 
that,  despited  the  disregard  of  nongeostrophic  advection  and  other 
important  physical  processes,  quasi-geostrophic  analysis  is  able  to 
resolve  features  having  a  much  smaller  scale  than  that  appropriate  to 
the  theory.  It  will  also  compare  vertical  motion  fields  calculated  by 
the  kinematic  and  quasi-geostrophic  methods,  and  will  attempt  to  relate 
fields  of  low  Richardson  number  to  the  clouds  and  circulation. 

The  primary  purpose  of  this  research  is  to  verify  and  refine 
the  model  of  the  three-dimension  flow  near  developing  comma  clouds  on 
the  basis  of  fine-scale  data,  so  that  forecasters  may  eventually  judge 
more  accurately  the  location  and  strength  of  disturbances  indicated  by 
this  type  of  characteristic  cloud  pattern. 
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CHAPTER  II 


CASE  SELECTION  AND  METHODOLOGY 


This  thesis  will  use  international  system  units  throughout 
except  for  pressure  (millibars),  vertical  motion  (mlcrobars/second),  and 
altimeter  setting  (Inches  of  mercury). 

A.  Case  Selection 

Inspection  of  geostationary  satellite  imagery  reveals  that  in 
only  one  case  during  the  1978  and  1979  Atmospheric  Variability 
Experiment  (AVE)  did  a  clearly  outlined  comma  cloud  cross  a  dense  rawin- 
sonde  network,  on  10-11  April  1979.  The  "regional  scale"  rawinsonde 
network  captured  all  of  the  associated  storm  at  lower  level s — the  dry 
southwest  quadrant,  the  moist  southeast  quadrant,  and  the  cool  airmass 
north  of  the  warm  front.  In  mid  to  high-levels,  the  trough  was  slightly 
west  of  the  AVE-SESAME  1  rawinsonde  network  until  0600  GMT  11  April  1979, 
but  standard  upper  air  data  supplemented  the  analysis  every  twelve  hours. 


Associated  tornadoes,  duststorms,  and  other  severe  weather 


occurred  In  the  SESAME  network  on  this  day.  Including  the  Infamous 
Wichita  Falls,  Texas,  tornado. 

Normally,  the  term  "comma  cloud"  refers  to  a  synoptic-scale 
cloud  shield  which  evolves  rather  gradually  In  a  developing  extratropi- 
cal  cyclonic  disturbance,  or  to  a  meso-ac  scale  cloud  In  the  cold  core  of 
an  upper  cyclone.  This  comma  cloud  was  unusual  In  that  It  developed 
very  rapidly  from  meso-c<  to  macro  proportions  above  an  intense  squall 
line  (in  the  middle  of  the  SESAME  data  set).  This  case  Is  represen¬ 
tative,  in  the  vigorous  extreme,  of  comma  clouds  which  accompany  rapid 
development  and  changing  weather,  and  a  detailed  knowledge  of  Its  struc¬ 
ture  is  valuable  to  operational  forecasters. 

B.  Data  Preparation 

Data  description.  Rawinsonde  data  were  available  every  three 
hours  from  1200  GMT  10  April  1979  to  1200  GMT  11  April  1979,  over  a 
"regional  scale"  network  with  approximately  214  km  average  separation  of 
each  station  from  Its  nearest  neighbor.  Service  C  teletype  rolls 
yielded  additional  12-hourly  reports,  while  Service  A  rolls  contained 
hourly  surface  observations.  Figure  1  Illustrates  the  SESAME  stations, 
the  standard  rawinsonde  stations  used,  and  the  analysis  grid  chosen. 

Half-hourly  satellite  photographs  consisted  of  half-mile 
resolution  video  data  In  the  daytime  and  one  mile  resolution  Infrared 
data  at  night,  with  the  National  Earth  Satellite  Service  MB  enhancement 
of  the  gray  scale.  Manually  digitized  radar  data  were  used  to  locate 


convective  features  when  grldding  was  missing  or  Incorrect  on  the  sat¬ 
ellite  images.  Nephanalyses  were  prepared  by  hand  from  these  satellite 
Images  and  radar  charts,  emphasizing  convective  areas,  the  duststorm, 
and  the  thunderstorm  cirrus  canopy  which  evolved  into  a  comma  cloud. 

Subjective  error  analysis.  Before  archiving  the  information 
on  magnetic  tape,  NASA  contractors  at  Texas  A&M  University  decoded  it, 
checked  it  for  errors  both  manually  and  by  machine,  and  interpolated  it 
linearly  from  the  surface  to  25mb  at  multiples  of  25  millibars 
(Fuelberg  (1974),  Gerhard  et  al .  (1979)).  Based  on  subjectively  ana¬ 
lyzed  constant-pressure  height  fields  for  each  time,  Gerhard  et  al . 
(1979)  decided  that  some  specific  rawlnsondes  were  probably  In  error 
because  their  heights  did  not  appear  to  fit  the  overall  height  pattern 
at  that  particular  time  (Sanders  (1981b)).  This  approach  to  error  iden¬ 
tification  recognizes  the  potential  for  error  In  every  Individual  rawln- 
sonde  ascent,  and  It  smooths  mesoscale  features  whether  they  represent 
noise  or  the  true  state  of  the  atmosphere.  However,  one  may  question 
whether  any  analyst  can  faithfully  distinguish  noise  from  real  data 
during  such  a  disturbed  weather  situation.  Moreover,  for  this  research, 
information  about  features  on  scales  down  to  480  kilometers,  persistent 
on  the  order  of  10  hours,  was  considered  acceptable.  It  was  more 
feasible  to  eliminate  any  errors  consistent  In  all  data  from  a  given 
station  (e.g.,  heights  reported  at  Gage,  Oklahoma  (GAG)  were  too  large 
at  all  times). 

For  this  reason,  twelve-hour  time-averaged  constant-pressure 
height  charts  were  prepared  for  700mb,  500mb,  and  200mb.  The  average- 
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height  charts,  subjectively  analyzed  at  a  20  meter  Interval,  are  shown 
as  Figures  2-4.  These  fields  suggested  that  several  stations.  Including 
Gage,  had  height  errors  persistent  in  time  and  pressure  level.  Although 
many  Interpretations  of  Figures  2-4  are  possible,  the  desired  result  of 
any  data  adjustment  must  be  to  weaken  the  height  gradients  and  curvature 
of  contours  which  the  original  data  Indicate  are  too  Intense.  This  was 
done  by  subtracting  a  constant  from  the  heights  of  all  rawlnsonde 
reports  for  the  three  most  erroneous  stations  (Gage,  Oklahoma  (17m 
high);  Abilene,  Texas  (l&n  high);  and  Stephenville,  Texas  (7m  high)), 
and  by  adding  a  constant  to  the  heights  of  a  fourth  station  (Goodland, 
Kansas  (5m  low)),  before  input  to  the  objective  analysis. 

Errors  In  station  pressure  can  lead  to  errors  in  the 
hydrostatically  Integrated  heights  above  the  surface,  and  Barnes  (1981) 
has  cast  doubt  on  the  dependability  of  station  pressure  for  the  special 
SESAME  rawlnsonde  stations.  However,  Barnes'  surface  pressure  analyses 
tend  to  remove  transient  features  in  the  same  manner  as  the  height  ana¬ 
lyses  by  NASA  and  the  persistent  errors  one  might  expect  In  the  height 
field  aloft  can  be  assessed  better  right  there--1n  the  height  field 
aloft  (Sanders  (1981b)).  Since  this  is  the  intent  and  effect  of  the 
corrections  to  heights  at  four  stations,  discussed  above,  Barnes'  sta¬ 
tion  pressure  error  estimates  were  not  used.  The  identification  and 
correction  of  Individual  errors  Is  explained  In  the  next  section  of  this 
chapter. 

Balloon  drift  and  release  time.  The  SESAME  data  sets  Include 
the  balloon  release  time  and  location  of  the  balloon  relative  to  the 


tracking  station  during  ascent.  Some  researchers  (e.g.,  Moore  &  Fuelberg 
(1981),  Jedlovec  &  Fuelberg  (1981))  have  used  this  Information  to 
Increase  the  accuracy  of  their  analyses.  In  this  study,  the  actual 
balloon  location  downwind  of  the  rawlnsonde  station  was  used,  rather 
than  station  location.  All  data  were  linearly  Interpolated  or  extrapo¬ 
lated  to  a  common  balloon  release  time  as  explained  In  Appendix  A. 

Static  stability.  The  static  stability  parameter  was  calcu¬ 
lated  at  each  station  at  each  time,  every  lOOmb  from  lOOOmb  to  lOOmb, 
with  the  following  equation: 

287 xT  x  { 1rt(e2>  -  ln(0i)}  (II.l) 

0  *  ■ 

p  x  50  mb 

where  T  is  the  observed  temperature  at  pressure  level  p,  and  83  and 
01  are  the  potential  temperatures  25mb  respectively  above  and  below 
pressure  level  p.  The  area  averaged  sigma  is  shown  in  Table  1  as  a 
function  of  time  and  pressure. 


Tlme(GMT) 


p(mb) 

NACA 

1200 

1500 

1800 

2100 

0000 

0300 

0600 

0900 

1200 

100 

1.830 

1.605 

1.619 

1.586 

1.628  1.617 

1.633 

1.720 

1.645 

1.639 

200 

.450 

.348 

.384 

.392 

.378 

.362 

.348 

.343 

.318 

.320 

300 

.070 

.060 

.061 

.070 

.062 

.061 

.058 

.057 

.059 

.077 

400 

.042 

.023 

.023 

.022 

.025 

.026 

.022 

.026 

.  .026 

.027 

500 

.027 

.021 

.023 

.019 

.019 

.021 

.017 

.019 

.019 

.019 

600 

.020 

.017 

.016 

.017 

.016 

.013 

.013 

.013 

.014 

.016 

700 

.015 

.017 

.014 

.016 

.017 

.016 

.014 

.014 

.017 

.014 

800 

.012 

.021 

.021 

.017 

.019 

.017 

.017 

.019 

.019 

.018 

900 

.010 

.035 

.037 

.035 

.021 

.021 

.024 

.025 

.021 

.023 

Table  1:  Area  averaged  static  stability  parameter  r  in  units  of  (mtykg) 
per  lOOmb  for  NACA  standard  and  for  every  three  hourly  SESAME  data  set. 

The  stability  parameters  for  NACA  Standard  Atmosphere 
(Haltlner  and  Martin  (1957))  were  computed  from  data  at  25mb  above  and 


below  each  pressure  level  of  Interest.  The  strength  of  the  maximum  of 
stability  at  900mb  Is  worthy  of  note. 


Auxiliary  soundings.  Because  of  the  location  of  the  cloud 
vortex  on  the  western  edge  of  the  SESAME  network,  it  was  necessary  to 
use  standard  12-hourly  rawlnsonde  data  from  stations  near  the  western 
part  of  the  grid  (see  Figure  1).  One  must  remember  that  the  analyses  at 
1200  GMT  10  April  1979,  and  0000  and  1200  GMT  11  April  1979  had  the 
benefit  of  this  extra  Information,  while  the  analyses  at  1500,  1800, 


2100,  0300,  0600,  and  0900  GMT  did  not. 


Height,  wind  direction,  and  wind  speed  for  these  12-hourly 
rawlnsondes,  were  extracted  from  Service  C  teletype  rolls  and  inter¬ 
polated  to  non-mandatory  levels  (1000,  900,  800,  and  600  mb).  A 
logarithmic  Interpolation  was  used  (see  Appendix  A  for  details). 

Richardson  number.  A  gradient  Richardson  number,  Rf,  was 
calculated  over  25mb  Intervals,  every  25mb,  for  each  rawlnsonde  by  the 
formula: 


R1 


9  a©  atu!  2 

0  az  9  Z 


(II. 2) 


where  g  Is  the  acceleration  of  gravity,  91s  the  potential  tem¬ 
perature,  |  u|  Is  the  wind  speed,  and  Z  Is  the  height.  Wind  speeds  and 
potential  temperatures  were  not  smoothed  or  averaged  in  the  vertical. 
Rawlnsondes  with  low  R^  (0<R^.25  and/or  0<Rj$.95)  In  two  vertical 
domains  (at/above  500mb  and  below  50Qnb)  were  plotted  on  charts  for 
every  three  hour  data  set.  Two  sets  of  charts  (one  above  500mb,  one 
below)  were  analyzed  by  hand  for  areas  where  0<R^.25  enclosed  In  areas 
where  0<Rj<.95. 


The  Richardson  number  gives  a  measure  of  the  ratio  of  buoyant 


forces  to  wind  accelerations  In  a  vertical  column  of  the  atmosphere. 
As  the  column  becomes  more  buoyant  ()6/aZ  decreases)  or  the  vertical 
difference  In  kinetic  energy  ((t|ui/*Z  )*)  Increases  Rj  decreases. 
Then  the  atmosphere  Is  more  likely  to  overturn  In  turbulent  exchange. 

Stone  (1966)  demonstrated  that  three  basically  different 
types  of  Instability  may  dominate  a  barocllnlc  current  under  different 
combinations  of  zonal  (k)  and  meridional  (X)  wavenumber,  and  of 
Richardson's  number  (R<).  The  problem  of  the  existence  and  relative 
power  of  each  of  these  forms  of  Instability  for  arbitrary  combinations 
of  values  of  K,  \  ,  R^  Is  far  from  solved,  but  Stone's  partial  solution 
suggests  very  strongly  that  one  may  expect: 

1.  Conventional  barocllnlc  instability  for  Ri>0.95. 

2.  Symmetric  Instability  for  0.25<Rj<0.95,  exciting  motion 
around  axes  parallel  to  the  zonal  flow. 

3.  Kel vln-Helmholz  instability  for  0.25<R{,  causing  motion 
around  horizontal  axes  perpendicular  to  zonal  flow. 

C.  Objective  Analysis 

Overview.  Cressman's  (1959)  successive  approximation  proce¬ 
dure  for  Interpolating  values  of  randomly  spaced  data  to  regularly 
spaced  gridpoints  was  used  In  this  study.  This  procedure  Is  sensitive 
to  many  factors  which  can  affect  the  quality  of  the  analysis,  such  as 
the  choice  of  the  rectangular  grid  spacing,  d,  and  the  maximum  radius, 
R,  at  which  data  points  Influence  a  grldpolnt.  An  overview  of  the 
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problems  Involved  will  make  the  solutions  taken  In  this  research  more 
understandable. 

Given  data  separated  In  the  average  by  s=214km,  one  would  like 
to  define  grldded  values  at  an  Interval  much  smaller  than  s.  However, 
the  expected  variance  of  Interpolation  error  (for  the  best  analysis 
conditions)  Increases  rapidly  as  the  grid  distance  decreases  below  s 
(Stephens  and  Stitt  (1970),  figure  9).  The  choice  of  a  grid  Interval 
only  slightly  less  than  s  ( d*s/l .33)  avoids  this  problem,  so  d  was  set 
to  161  km.  Stephens  and  Stitt  (1970)  also  Identified  a  range  of  radii 
of  Influence,  R  (defined  below)  for  Cressman  analysis  of  data  with 
separation  s  to  a  grid  of  Interval  d,  over  which  the  error  generated  by 
the  Interpolation  operator  Is  minimal.  The  strategy  of  using  suc¬ 
cessively  smaller  R  with  Increasingly  more  restrictive  error  iden¬ 
tification,  must  be  done  within  this  range.  Stephens  and  Stitt  advise 
avoiding  the  lower  limit  on  R,  approximately  1.5d,  as  deduced  from  their 
figure  9.  Preliminary  runs  with  values  of  R  from  1.5d  to3.5d  Indicated 
that  this  SESAME  data  set  yields  a  reasonably  smooth  height  analysis  for 
R>3d. 

Table  2  (below)  shows  the  five  steps  which  were  applied  In 
this  analysis,  with  R  ranging  from  3d  down  to  2d.  Pass  I  created  a 
first  guess  field  from  the  reported  data,  without  reference  to  data  at 
any  other  pressure  level  or  time.  Every  field  except  the  last  pass  was 
smoothed  with  a  simple  smoothing  operator  (see  Appendix  B)  to  eliminate 
high  frequency  noise,  although  the  Cressman  weighting  function  filters 
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much  noise  at  wavelengths  less  than  2R  (Stephens  1967)).  No  analysis, 
objective  or  subjective,  can  extract  Information  about  disturbances  of 
wavelengths  less  than  2s,  the  lower  limit  of  the  ability  of  the  SESAME 
data  set  to  discriminate  patterns  unambiguously.  Therefore,  for  the 
last  pass,  when  2R=3s  (since  R=2d  and  d=s/1.33),  noise  at  wavelengths 
less  than  the  2s  lower  limit  of  the  SESAME  data  set  was  definitely  eli¬ 
minated  by  the  Cressman  analysis  without  need  for  final  smoothing. 

Since  the  observations  may  contain  unacceptlble  errors,  the  Pass  I  field 
was  used  only  to  Identify  the  most  erroneous  data  (by  a  method  which 
will  be  discussed  below).  After  Passes  II,  III,  and  IV,  the  latest 
values  of  all  fields  were  substituted  for  missing  reports  in  the  origi¬ 
nal  data  set  and  used  as  bogus  for  the  next  pass. 


Pass  I 

R*3.00d  *  483km 

Smooth. 

Check  errors,  but 

do  not  bogus. 

Phase 

II 

R*3.00d  -  483km 

Smooth. 

Check  errors  and  bogus  missing  data. 

Phase 

III 

R-2.67d  *  429km 

Smooth. 

Check  errors  and  bogus. 

Phase 

IV 

R=2.33d  =  376km 

Smooth. 

Check  errors  and  bogus. 

Phase 

V 

R*2.00d  *  322km 

Do  not  smooth. 

Table  2:  Analysis  strategy. 

Computational  grid.  The  analysis  region  has  11x12  points  In 
the  horizontal  and  eleven  levels  In  the  vertical:  the  surface  plus 
every  lOOmb  from  lOOOmb  to  lOOmb.  Gridpolnts  below  the  smooth  terrain 
analysis  were  flagged  Immediately  after  the  height  analysis  for  each 
data  set.  The  pressure-coordinate  scheme  Is  simpler  than  a  terrain 
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following  ^-coordinate  system  because  the  data  were  supplied  at 
constant  pressure  levels.  The  system  would  require  Interpolation  of 
data  from  pressure  to  ^-coordinates  for  computation  then  back  to 
pressure  coordinates  for  graphic  display.  It  Is  hoped  that  vertical 
finite  differencing  over  Intervals  of  constant  pressure  would  yield 
more  consistent  results  than  differencing  between  layers,  which  vary  In 
pressure  separation.  For  these  reasons,  the  pressure  coordinate  system 
was  chosen. 

Cressman  analysis.  Data  within  a  fixed  horizontal  "radius  of 
Influence",  R,  of  a  given  grldpolnt  were  averaged  after  weighting  by 


the  Cressman  (1959)  weighting  function,  w: 

n  n 

So  *  Z  ("ISO  /Z  *1 

1=1  1-1 


for  wi’ 


R2-di2 


(II. 3) 


(II.4) 


R2+di2 

where  data,  ,  observed  at  locations  (x^,  y^)  within  the  radius  of 
Influence,  R,  of  the  grldpolnt  (x0,  y0)  were  weighted  relative  to  the 
distance  df.  Each  pressure  level  was  considered  independent  of  data 
above  or  below  the  level.  If  rawlnsonde  reports  of  height  and  wind 
were  both  available,  the  height  estimate,  Zj ,  was  augmented  by  the 
height  gradient  computed  from  the  observed  wind  assuming  geostrophlc 
balance,  and  the  compound  report  was  weighted  four  times  a  simple 
height  report  (Cressman  (1959)).  Appendix  B  details  further  theory 
behind  this  scheme. 

Objective  error  Identification.  After  each  analysis  pass, 
all  data  were  checked  against  an  analyzed  value  at  the  rawlnsonde  loca- 


tlon,  calculated  by  bilinear  interpolation  of  the  latest  gridded  analy¬ 
sis  (see  Appendix  B).  Each  deviation  was  checked  against  an  error  cri¬ 
terion  which  varied  with  pressure  level,  with  parameter,  and  with  the 
radius  of  Influence,  R.  If  the  deviation  was  less  than  the  criterion, 
the  observed  value  was  stored  in  a  working  rawinsonde  array. 

Otherwise,  the  observation  was  disregarded  for  that  particular  pass, 
and  the  analyzed  value  was  stored  in  the  working  array.  This  working 
rawinsonde  array  was  used  to  calculate  the  next  pass  analysis.  The 
original  reports,  including  missing  flags,  were  stored  for  reference 
throughout  the  analysis,  and  re-checked  against  later  pass  analyses. 
Thus,  observations  considered  "bad"  relative  to  an  early,  rough  analy¬ 
sis  could  be  used  later  when  their  comparison  to  more  reasonable  analy¬ 
ses  indicated  the  observation  was  "good." 

A  visual  Inspection  of  the  observation/analysis  deviations, 

|  z0-za|  and  ((uo-ua)2+(vo-va)2)1^.  f°r  some  of  the  data  showed  that  the 

statistical  properties  of  the  data  set  itself  should  dictate  the  height 

and  wind  error  criteria.  Each  of  the  nine  data  sets  were  analyzed  by  a 

single  pass  with  R=3d,  giving  nine  separate  analyses  of  z,  u,  and  v  at 

each  level,  and  the  root-mean-square  (r.m.s.)  of  all  data  at  each 

pressure  level  was  accumulated  and  are  illustrated  in  Table  3. 

p(mb)  sfc  1000  900  800  700  600  500  400  300  200  100 

oz(m)  -  8.5  14.8  16.9  15.1  18.8  19.6  22.6  28.5  30.5  24.7 

ov(m/s)  4.8  3.6  5.4  6.6  6.5  7.3  8.8  9.8  11.3  9.2  6.5 

Table  3:  Standard  deviation  of  first-pass  analyzed  heights  and  wind 

speeds  accumulated  from  all  nine  sets  of  observed  data. 
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These  standard  deviations  were  used  to  define  basic  height- 
dependent  criteria, $z  and  $v,  for  testing  errors  in  height  and  wind, 
respectively.  The  basic  criteria,  and  £r,  for  temperature  and  rela¬ 
tive  humidity  were  taken  as  twice  the  r.m.s.  errors  expected  for  AVE 
data  (Fuelberg  (1974)),  independent  of  height.  For  the  surface 
pressure  tendency,  which  is  used  as  part  of  the  lower  boundary  con¬ 
dition  on  kinanatlc  omega,  8%  is  taken  as  50X  in  order  to  eliminate 
only  the  worst  errors. 

In  order  to  restrict  the  error  tolerances  more  and  more  with 
each  analysis  pass,  the  criteria  for  error  Identification  were 
decreased  linearly  in  R  from  2ST(p)  at  R=3d,  to  ^(p)  at  R=2d,  where 
8r  is  one  of  the  basic  pressure-dependent  criteria  relative  to  the 
parameters  (ir)  which  were  analyzed.  Table  4  indicates  the  numerical 
definition  of  the  basic  criteria  and  the  testing  inequalities  which 
were  applied  to  each  piece  of  data  after  each  analysis  pass. 

Subscripts  o  and  a  refer  to  observed  data  and  analyzed  values  of 
height,  u  and  v  wind  components,  pressure  tendency,  temperature  and 
relative  humidity;  <TZ  and  are  pressure-dependent  standard 
deviations  of  height  and  wind  speed  for  the  entire  SESAME  data  set. 

The  factor  (R-l)  decreases  the  criteria  linearly  with  R,  the  radius  of 
Influence. 


I 


II 


1 

i 


z 

(height) 


Basic  Criterion: 
6z=oz(p)(meters) 


Application: 

lZ0-Zal>(R-l)6z 


^  6v=Oy{p)(meters/seCOnd) 

(vector  wind) 


l(u0-ua)2+(vo-va>2l1/2>(R-l>6, 


x  6x=50% 

(surface  pressure  tendency) 


(xo‘Xa^x0MR-U6x 
or  (xoxa)<0 


T  6t=2°C 

(temperature) 


IVTal>(R-l)«i 

|rh0-rha|/rh0  >(R-l)<5r 

Table  4:  Criteria  (5Z,  &v»  ^T»  ^r)  for  error  identification.  See 

text  for  explanation.  Also  see  Fuelberg  (1974)  for  and  £r. 


rh  6r=20% 

(relative  humidity) 


h 


i 


i 


Residual  errors.  Table  5,  below  shows  the  r.m.s.  deviations 
of  analyzed  fields  from  the  observed  data  after  five  analysis  passes, 
with  a  final  radius  of  influence  R=2d.  Values  at  other  pressure  levels 
are  listed  in  Tables  7  and  8  in  Appendix  B. 


Time  (GMT) 

z(m) 

|t| (m/s) 

z(m) 

It  1 (m/s) 

(500mb) 

(200mb) 

1200 

5.5 

3.1 

9.6 

5.3 

1500 

6.6 

3.2 

9.8 

6.3 

1800 

5.5 

3.4 

12.1 

4.5 

2100 

6.5 

5.0 

39.8 

5.3 

0000 

8.8 

4.9 

23.8 

3.8 

0300 

6.9 

4.8 

23.1 

4.8 

0600 

6.4 

8.2 

16.6 

9.1 

0900 

10.0 

6.4 

14.0 

13.2 

1200 

11.4 

4.0 

10.5 

6.5 

(S2=19.6) 

{ 8 .8*6  y ) 

(6Z=30.5) 

(9.2=6 v) 

Table  5:  Final  r.m.s.  height  and  vector  wind  speed  errors  for  each 


data  set  at  500mb  and  200mb. 


ti 
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Some  of  the  values  at  2O0mb  are  larger  than  the  criteria 
5Z  and  Sv.  and  they  deserve  the  following  explanation. 

At  2100  GMT,  <Ti(200mb)=18.6m  for  31  stations  inside  the  ana¬ 
lysis  grid.  However,  one  station  (COU)  with  a  200m  indicated  error 
raises  the  Tz  for  all  32  stations  to  39.8m.  The  analyzed  height  at  COU 
is  consistent  in  time  and  space  with  all  data,  except  the  2100  GMT 
sounding  at  COU.  The  COU  heights  are  reasonable  as  high  as  500  mb,  but 
the  indicated  error  increases  from  38m  at  400mb,  to  115m  at  300mb,  to 
200m  at  200mb;  and  they  lose  consistency  with  the  1800  and  0000  GMT 
soundings  from  COU.  Temperature  reports  are  exceedingly  warm  in  this 
part  of  the  sounding.  All  the  evidence  indicates  that  the  2100  GMT 
rawinsonde  from  COU  failed  above  500mb,  that  the  reported  heights  were 
bad,  and  that  the  analysis  scheme  was  able  to  identify  the  error. 

AT  0900  GMT,  rv(200mb)=13.2,  mainly  because  the  analysis 
smoothed  a  very  strong  jet  maximum.  The  strongest  wind  at  200mb  at  that 
time  was  at  OKC,  reported  at  83  m/s  but  analyzed  as  41  m/s,  with  an  43 
m/s  vector  deviation  from  the  observation. 

A  thorough  search  for  the  causes  of  large  r.m.s.  errors  indi¬ 
cated  that  the  height  analysis  dealt  with  height  errors  very  well,  but 
that  the  wind  analysis  smoothed  some  jet  streaks  to  half  the  intensity 
a  hand  analysis  shows  and  also  removed  mesoscale  troughs  evident  in 
hand-drawn  500mb  streamline  analyses.  From  the  understatement  of  hori¬ 
zontal  wind  gradients  in  such  regions,  one  must  expect  weaker  values  of 
kinematically  computed  vertical  motion  than  a  better  wind  analysis 
might  reveal.  The  results  from  this  kinematic  analysis  will  be  com- 
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pared  to  kinematic  omega  computed  by  Moore  and  Fuel  berg  (1981)  for 
reasonableness. 


0.  Vertical  Velocity 

Comparison  of  methods.  Two  methods  of  calculating  vertical 
motion  were  used  in  this  study,  namely  the  kinematic  and  the  quasi- 
geostrophic  (QG)  methods.  The  QG  omega  equation  was  chosen  for  the  two 
following  reasons.  First,  the  cause  of  vertical  motion  may  be 
separated  into  two  physically  meaningful  parts.  The  QG  analysis  is  a 
compromise  between  the  kinematic  method,  which  provides  no  such  physi¬ 
cal  insight,  and  much  more  complicated  omega  equations  which  also 
account  for  latent  heat  release,  surface  friction,  advection  by  the 
divergent  part  of  the  wind,  and  so  on.  Second,  these  calculations, 
based  entirely  on  a  smooth  height  analysis,  eliminate  the  large  cumula¬ 
tive  error  characteristic  of  the  kinematic  method  which  can  result  from 
minor  noise  in  the  wind  field.  Both  of  these  points  present  paradoxes, 
however,  which  underscore  the  independence  of  the  two  methods. 
Krishnamurti  (1968)  demonstrated  what  operational  forecasters  have  long 
known,  that  the  two  first-order  terms  of  QG  forcing  do  not  always  ade¬ 
quately  define  the  total  vertical  motion  of  the  atmosphere  in  response 
to  complex  forcing.  Using  only  the  first  two  terms,  one  must  rely  on 
experience  and  empirical  guidance  to  decide  whether  the  non-QG  forces 
may  dominate  in  some  area,  and  what  changes  they  will  cause  in  the 
height  field.  Second,  the  kinematic  method  (adjusted)  often  yields 
omega  fields  which  compare  well  to  the  observed  clouds  and  weather 
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despite  the  problem  of  noise.  Again,  the  analyst  can  never  be  certain 
whether  any  particular  kinematic  analysis  Is  better  or  worse  than  the 
corresponding  QG  analysis  without  comparing  the  fields  to  cloud,  preci¬ 
pitation,  or  moisture  patterns,  or  by  calculating  some  higher  order 
version  of  the  omega  equation  (e.g.,  Krishnamurti  (1968)). 

Summary  of  computations.  With  the  quasi geostrophic  (QG) 
omega  equation,  vertical  motion  was  calculated  over  a  7x8  gridplnt 
Interior  domain  (<u>=o  around  around  the  sides)  using  boundary  conditions 
specified  at  lOflmb  (c*N>)  and  at  the  first  vertical  level  above  terrain 
(generally  900mb  In  the  east  and  80Qnb  in  the  west).  Kinematic  omega 
was  computed  over  a  9x11  grid  and  adjusted  to  zero  at  lOOmb  using 
O'Brien's  (1969)  quadratic  technique.  Two  different  types  of  lower 
boundary  condition  were  used:  one  for  the  kinematic  method  and  another 
for  the  QG  method. 

The  QG  omega  field  was  partitioned  in  two  separate  ways.  The 
traditional  partition  between  forcing  due  to  differential  vorticity 
advection  and  to  thermal  advection,  and  also  an  alternate  partition 
between  forcing  due  to  advection  of  vorticity  by  the  thermal  wind  and 
to  a  deformation  function  both  were  computed. 

Top  boundary  conditions.  At  the  top  (lOOmb)  omega  can  be 
assumed  to  be  zero  or  it  can  be  calculated  by  some  Independent  method, 
e.g.  the  adiabatic  method.  Omega  was  calculated  at  lOOmb  by  the  (dry) 
adiabatic  method  In  order  to  Inspect  the  magnitude  of  stratospheric 
vertical  motion. 
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Calculations  by  this  method  indicated  values  of  omega  at 
lOOmb  in  the  range  +1.15^*b/s  for  all  of  the  nine  data  sets  (omega 
fields  not  illustrated).  Thus  a  representative  value  of +^.5^*b/s  would 
be  susceptible  to  around  +^17 ^*b/s  rms  error,  as  discussed  in  Appendix 
C.  As  might  be  expected,  this  magnitude  of  error  resulted  in  erratic 
continuity  in  the  computed  omega  fields.  Therefore,  the  top  boundary 
condition  was  assumed  to  beco=o» 

Lower  boundary  conditions.  Two  types  of  lower  boundary  con¬ 
ditions  were  used,  one  for  each  method.  First,  the  kinematic  method 


requires  a  surface  vertical  motion  which  results  from  observed  wind: 


“s*  "  Ps9Ws#7h)  +  xs 


(II. 5) 


where  j>s,  Vs,  and  ”XS  are  the  density,  observed  wind,  and  local 
pressure  tendency  at  the  surface,  h  is  the  terrain  height,  and  g  is  the 
acceleration  of  gravity.  Surface  pressure,  density,  and  geostrophic 
wind  (used  in  another  boundary  condition)  were  computed  indirectly 
based  on  the  hypsometric  equation,  as  detailed  in  Appendix  C.  The 
observed  wind  V^s  was  gadded  from  rawinsonde  u  and  v  components, 
without  the  use  of  hourly  surface  reports. 

The  QG  omega  equation  is  most  consistent  with  the  lower  boun¬ 


dary  condition: 


Here,  the  surface  geostrophlc  wind,  V^g,  Is  taken  as  half  the  value 
Indicated  by  the  gradient  of  pressure,  P,  at  the  terrain  height  of  each 
grldpolnt.  In  order  to  account  for  the  decrease  In  wind  speed  toward 
the  surface  because  of  friction  (Carr  (1981)).  See  Appendix  C  for  the 
method  of  computing  P  hydrostatically  from  the  gridded  height  field. 

Kinematic  omega.  The  kinematic  method  for  calculating  omega 
consists  of  integrating  the  continuity  equation  from  the  surface 
upward. 

N  i 

“N*“o  +  l  (pk-l-pk)  5  («k-l+«k)  (II,8) 

k=l  c 

The  average  divergence,  ^(Sk-l+Sjj) t  over  a  pressure  layer,  p^.i  to  p^, 
determines  the  vertical  change  In  omega  over  the  layer.  Only  a  lower 
boundary  value,  t»fc,  Is  required. 

If  one  assumes  wind  errors  at  neighboring  gridpolnts  to  be 
independent  of  one  another,  then  the  horizontal  gradient  of  wind  would 
have  an  r.m.s.  error  of  <Tv/d,  because  the  winds  at  either  end  of  the  2d 
interval  would  each  contribute  o'*  error.  As  shown  in  Appendix  C,  this 
leads  one  to  expect  errors  of  about  10 ^»b/s  in  omega  after  vertical 
integration  of  divergence.  The  r.m.s.  value  of  omega  at  lOOmb  calcu¬ 
lated  by  the  kinematic  method  before  correction  to  zero  was  slightly 
smaller  than  this  expectation  (~7.5^*b/s). 

A  simple  but  effective  technique  for  adjusting  the  solution 
of  the  kinematic  equation  (II. 8)  to  some  reasonable  top  boundary  value 
Is  that  due  to  O'Brien  (1969).  The  top  value,  ft)*,  of  a  column  of  K+l 
calculated  values  of  omega  Is  In  error  by  an  amount  (w^- wtop)  from  the 
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desired  boundary  value,  wtAn*  One  can  distribute  a  correction  to  each 

k*0,K  (II. 9) 

k*0,K  (II. 10) 

Experimentation  Indicated  that  equation  {II. 9)  distributes  too  much  error 
into  the  middle  part  of  the  column,  creating  bulls-eyes  In  the  omega  field 
in  mid-levels,  more  often  than  not  of  the  opposite  sign  one  would  expect 
from  synoptic  reasoning.  The  quadratic  adjustment  (II. 10)  controls  the 
corrections  at  lower  levels,  eliminates  the  bulls-eyes  at  500mb,  and 
gives  somewhat  satisfying  results  at  high  levels. 

In  summary,  the  final  product  of  O'Brien's  quadratic  tech¬ 
nique  in  this  research  conforms  very  well  to  the  synoptic  pattern, 
suggesting  that  the  actual  noise  problem  is  less  severe  than  one  might 
expect  from  the  theory  of  the  kinematic  method. 

The  quasi-geostrophic  omega  equation.  In  three  dimensions, 

V2u>  +  f02/o  =  F  (II. 11) 

ap2 

where  f0  is  a  constant  value  of  the  coriolis  parameter,  F(x,  y,  p)  is  a 
forcing  function,  and  static  stability,  <T  ,  depends  only  on  time  and 
pressure.  The  vertical  derivative  was  computed  using  finite  differences 
over  two  hundred  millibars. 


ttfct either  linearly  with  height: 
k  , 

“*k*“k  -  £  'uK"“top' 
or  quadratlcally: 

k(k+l) 

“  k  “k  -  (“K"utop' 


There  are  two  popular  ways  to  express  the  Q6  forcing  function: 


the  traditional 


F  -  -  - —  (-$g*U+f))+-  v*(-^.7T) 

o  3(-p)  op 

=  F1+F2 

and  the  alternative 

f0  f0  3D  3E 

F  =  -  (-VT*V(2c+f)]  -  —  (E—  -  D—  ) 

O  O  dp  dp 

=  F3+F4 


(11.12) 


(11.13) 


The  first  term  In  equation  (11.12),  FI,  represents  forcing  by  the  dif¬ 
ferential  advectlon  of  absolute  vorticity  (f+f)  by  the  geostrophic  wind 
$j.  The  second  term,  F2,  Is  forcing  by  the  laplacian  of  thermal  advec- 
tion  by  Vg.  The  alternative  partition  due  to  Winn-Nielsen  (1959) 
Involves  Trenberth's  (1977)  F3*(2A+C)  term  for  the  advection  of  the 
quantity  (2f  +f)  by  the  thermal  w1nd,7j,  and  Winn-Ni  el  sen's  F4=(-2A) 
deformation  function,  where  0  and  E: 


dv  du  d  u  dv 

0  *  —  ♦  — ;  E  - -  (11.14) 

dx  dy  dx  dy 


are  the  so-called  shearing  and  stretching  parts  of  deformation. 
Finite-difference  forms  and  their  mutual  consistency  are  discussed  in 
Appendix  C. 

Winn-Niel son's  formulation  permits  an  easy  visualization  of 
the  Q6  vertical  motion  from  an  isobarlc  chart  of  thickness  and  (2£+f), 
providing  the  F4  forcing  function  is  negligible  compared  to  F3.  Using 
the  thermal  wind  indicated  by  the  direction  and  spacing  of  thickness 
"isotherms"  (cold  on  the  left),  the  degree  of  cyclonic  vorticity  advec- 
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tlon  Is  Inversely  proportional  to  omega,  without  the  troublesome  * - 

i(-p) 

and  72  operators  In  the  FI  and  F2  terms.  One  of  the  secondary  objec¬ 
tives  of  this  research  Is  to  examine  the  magnitude  of  F4  and  the 
resulting  QG  omega  using  fine-scale  data.  In  order  to  verify  the  utility 
of  this  analysis  tool  In  the  vicinity  of  a  strong  developing  wave. 


CHAPTER  III 

RESULTS 

A.  Vertical  Motion 

A  few  words  of  explanation  will  facilitate  the  discussion  of 
the  nine  individual,  three-hourly  data  sets. 

The  clouds  and  weather  on  10  April  1979  were  influenced  very 
strongly  by  five  meso-<x  (see  Orlanski  (1975))  short  waves  that  passed 
through  the  area.  These  disturbances  had  wavelengths  from  500  to  800 
km,  phase  speeds  of  20  to  33m/s,  and  amplitudes  of  around  30m  in  the 
500mb  height  field.  They  can  be  followed  best  in  the  500mb  relative 
vorticity  analyses  which  varied  ±5  x  10"5s"l  as  the  minor  troughs  and 
ridges  passed.  Waves  of  this  scale  in  the  extratropical  easterlies  have 
not  received  much  attention  in  the  literature.  Because  these  short 
waves  are  smaller  than  the  better  understood  synoptic  scale  short  waves 
(Palmen  and  Newton  (1969))  they  will  be  called  "minor  short  waves"  or 
"minor  waves"  in  this  thesis.  Other  researchers  have  found  evidence  of 
these  disturbances  in  this  case.  They  show  up  in  the  300mb  isotach  ana¬ 
lyses  In  Moore  and  Fuelberg  (1981)  and  in  the  time  cross-sections  of 
weather,  humidity,  and  vertical  motion  in  Wilson  (1981). 


In  this  thesis,  500mb  ascent  downwind  of  subsidence  will  be 
taken  as  evidence  of  a  minor  wave  in  the  mid-  to  upper-troposphere  which 
can  be  corroborated,  hopefully,  in  the  height  and  vorticity  fields. 

Each  use  of  the  term  "500mb  jet"  in  this  discussion  will  refer 
to  the  zone  of  maximum  winds  implied  by  a  region  of  strong  500mb  height 
contour  gradient  in  Figures  5c,  6c,  7c,  etc. 

For  brevity,  the  total  quasi-geostrophic  omega  resulting  from 
the  sum  of  FI  and  F2  forcing,  will  be  called  Q6  omega  in  constrast  to 
the  four  partitions  of  this  field  by  differential  vorticity  advection 
(FI  omega),  laplaclan  of  thermal  advection  (F2  omega),  advection  of  vor¬ 
ticity  by  the  thermal  wind  (F3  omega),  and  Winn-Ni el  sen's  deformation 
function  (F4  omega).  "Kinematic  omega"  refers  to  omega  computed  kinema¬ 
tically  and  adjusted  by  O'Brien's  quadratic  method.  Omega  of  absolute 
value  greater  than  10y>/s  is  defined  as  "strong"  vertical  motion. 

Bear  in  mind  the  fact  that  these  data  sets  are  very  close  in 
time.  Features  in  the  omega  analyses  which  exhibit  time  continuity 
should  be  considered  true  expressions  of  the  atmosphere  as  interpreted 
by  the  respective  mathematical  model  (QG  or  kinematic),  but  features 
which  appear  and  disappear  within  six  hours  must  be  judged  as  noise. 

Reference  to  tornado  and  hail  occurrences  are  taken  from  the 
damage  summaries  in  Alberty,  et  al .  (1979). 

Daybreak,  10  April  1979.  At  1200  GMT  (Fig  5a  and  b)  showers 
accompany  a  low  In  the  lee  of  the  Colorado  Rockies,  a  cold  front  in  New 
Mexico,  and  a  developing  warm  front  on  the  Texas  coast.  Gulf  stratus 
extends  up  the  Rio  Grande  valley.  A  500mb  trough  in  Kansas  (Fig  5c)  is 
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minor  wave  #1  of  this  case  study.  A  500  mb  ridge  across  the  Texas 
panhandle  separates  this  trough  from  minor  wave  #2  in  New  Mexico. 

Both  the  QG  and  kinematic  omega  fields  (Fig  5d  and  e)  show 
ascent  In  western  Kansas  and  subsidence  in  Oklahoma  associated  with 
minor  wave  #1.  The  ascent  in  this  case  is  not  directly  downstream  from 
the  descent  but  Is  located  in  the  cold  air  on  the  left  of  the  500mb  jet 
which  flows  from  south  to  north  across  Kansas  (see  Fig  5c).  Kinematic 
omega  shows  500mb  ascent  in  the  showery  region  of  the  Texas  coast,  while 
QG  omega  does  not.  Both  methods  show  an  area  of  ascent  in  eastern  New 
Mexico  at  least  six  hours  before  the  appearance  of  a  subsynoptic  surface 
low  (SSL)  in  west  Texas  (see  Fig  8a). 

The  two  methods  disagree  about  the  situation  over  extreme  west 
Texas.  Kinematic  omega  (Fig  5e)  suggests  rising  motion  in  advance  of 
another  minor  wave,  and  kinematic  omega  three  hours  later  (Fig  6e)  re¬ 
inforces  this  Inference,  as  weak  subsidence  appears  south  of  El  Paso, 
Texas.  QG  omega  shows  moderate  subsidence  (+7^b/s  in  Fig  5d)  over 
southern  New  Mexico.  Surface  temperatures  at  Carlsbad  and  Roswell,  New 
Mexico,  rose  14°C  In  two  hours,  suggesting  warming  by  strong  subsidence. 
The  dry  line  probably  began  moving  some  time  In  the  next  three  hours, 
but  it  cannot  be  located  well  until  1500  GMT  (Fig  6a)  In  southeastern 
New  Mexico.  The  surface  warming  supports  the  QG  omega  analysis  of  sub¬ 
sidence  behind  minor  wave  #2. 

As  indicated  In  Fig  5f  and  g,  FI  and  F2  omega  generally  te 
to  cancel  one  another,  except  near  strong  centers  where  one  or  the  other 
often  dominates.  In  west  Texas  and  southern  New  Mexico,  cold  advection 
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Is  balancing  subsidence  warming,  so  that  500mb  temperatures  are  nearly 
constant  (Fig  5c  and  6c).  Although  one  might  expect  the  subsidence 
warming  to  weaken  faster  than  the  horizontal  advectlon  cooling  well 
below  500mb  (because  the  ground  does  not  interrupt  horizontal  motion  as 
directly  as  vertical),  with  consequent  low  level  cooling,  the  opposite 
occurs.  That  is,  surface  temperatures  are  rising  at  the  rate  of  14°C  in 
two  hours. 

Figures  5h  and  5d  show  that  F3  omega  approximates  (total)  QG 
omega  very  well,  but  exaggerates  the  maximum  value  by  25%.  For  reasons 
discussed  in  Appendix  C  (see  equation  C.28),  deformation  forced  omega 
(F4)  was  calculated  by  an  operator  inconsistent  with  the  residual 
(F1+F2=F3),  so  that  the  sum  of  Fig  5h  and  i  may  be  slightly  different 
from  the  QG  omega  of  Fig  5d. 

Mid-Horning,  10  April  1979.  By  1500  GMT  (Fig  6a  and  b)  the 
cold  front  is  slowing  down  in  New  Mexico  upwind  of  a  region  of  pressure 
falls  in  southeastern  New  Mexico.  A  dry  line  is  apparent  in  this  area 
and  in  the  Pecos  Valley  of  Texas.  Visual  satellite  data  shows  Gulf 
Coast  stratus  as  far  northwest  as  Lubbock,  Texas  (dewpoint  10°C). 

Despite  the  gradual  clearing  of  clouds  in  east  Texas,  surface  winds, 
temperatures,  and  dew  points  (not  illustrated)  show  a  strengthening  warm 
front  in  the  area  (Fig  6a).  Now  three  minor  waves  are  evident  in  the 
500mb  relative  vorticity  pattern  (Fig  6c):  #1  on  the  Kansas/Nebraska 

border,  #2  in  northwest  Texas,  and  a  new  minor  wave  #3  in  extreme  west 
Texas.  Wave  #2  seems  to  be  weakening. 


27 


l 


4 


4 


4 


4 


The  QG  and  Kinematic  omegas  (Fig  6d  and  e)  both  reflect  minor 
wave  #3  In  west  Texas  with  a  small  region  of  ascent  directly  over  the 
northern  part  of  the  surface  dry  line,  downwind  of  descent  near  El  Paso, 
Texas,  in  the  area  between  the  dry  line  and  the  cold  front.  Kinematic 
omega  shows  rising  motion  in  northeast  Colorado  and  subsidence  in  north¬ 
western  Oklahoma  evidently  associated  with  wave  #1. 

Kinematic  omega  shows  ascent  in  south  Texas  near  a  thun¬ 
derstorm  area  on  the  coast.  This  area  moves  eastward  over  the  following 
six  hours  (Fig  6e,  7e,  and  8e).  Weak  subsidence  covers  northeast  Texas 
in  Fig  6e,  and  skies  are  clearing  there  (Fig  5b,  6b,  and  7b). 

The  QG  omega  pattern  (Fig  6d)  shows  very  little  structure  on 
the  anticyclonic  side  of  the  500mb  jet,  which  flows  from  southwest  to 
northeast  across  west  Texas  and  Oklahoma  (Fig  6c),  especially  in  com¬ 
parison  with  kinematic  omega  (Fig  6e). 

Figures  6f  and  g  reveal  that  the  ascent  probably  associated 
with  the  nascent  SSL  in  eastern  New  Mexico/west  Texas  is  due  entirely  to 
FI  omega.  The  subsidence  in  extreme  west  Texas  behind  the  dry  line  is 
forced  as  much  by  vorticity  advection  as  by  cold  advection,  in  contrast 
to  the  domination  by  F2  omega  three  hours  earlier  (Fig  5f  and  g). 

Again  F3  omega  (Fig  6h)  approximates  total  QG  omega  (Fig  6d) 
closely  with  some  exaggeration  of  the  values  in  closed  centers. 

Noon,  10  April  1979.  At  1800  GMT  (Fig  7a  and  b)  severe  thun¬ 
derstorms  break  out  from  the  bulge  on  the  dry  line  northward  toward 
Amarillo,  Texas.  Severe  storms  break  out  also  in  north  central  Texas 
between  the  warm  front  and  the  Red  River.'  Wink,  Texas,  and  Carlsbad, 
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New  Mexico,  first  report  blowing  dust  at  1800  GMT.  Minor  wave  #2  is 
barely  evident  in  the  500mb  relative  vorticity  pattern  (Fig  7c)  in 
northwest  Oklahoma.  Minor  wave  #3  is  in  southeastern  New  Mexico  just 
upstream  from  the  first  (weak)  tornado  of  the  day  (near  the  thunderstorm 
symbol  in  west  Texas  on  Fig  7b).  Low  clouds  are  dissipating  in 
southeast  Colorado  after  at  least  six  hours  of  downslope  flow. 

QG  omega  (Fig  7d)  shows  moderate  uplift  at  the  dry  line  bulge. 
The  ascent  region  of  kinematic  omega  in  eastern  New  Mexico  failed  to 
move  from  1500  to  1800  GMT  (Fig  6e  and  7e),  but  remained  along  the  cold 
front,  where  a  few  thunderstorms  are  already  occurring  (Fig  7b).  This 
is  a  good  example  of  the  ability  of  the  kinematic  method  to  pick  out 
small  regions  of  active  convection.  The  kinematic  omega  ascent  region 
over  the  Red  River  is  remarkably  timely--after  many  hours  of  rain  and 
drizzle,  this  small  area  breaks  out  in  hailstorms  just  before  the  time 
of  this  analysis  (see  Alberty,  et  al .  (1979)).  QG  omega  does  not  show 

moderate  ascent  over  this  thunderstorm  area  until  2100  GMT  (Tig  8d), 
three  hours  after  kinematic  omega  does.  Although  the  lower  tropospheric 
convergence  resulted  in  a  strong  upward  flux  of  mass,  moisture,  and  sen¬ 
sible  heat,  the  transport  probably  was  confined  to  a  few  violent  thun¬ 
derstorm  updrafts.  As  a  result,  the  600  to  400mb  height  patterns  and, 
hence,  the  500mb  QG  omega  pattern  remained  unchanged  on  the  scale  of 
this  analysis  for  the  first  few  hours  of  this  convection. 

In  east  Texas,  QG  omega  subsidence  disagrees  with  kinematic 


omega  rising  motion  (Fig  7d  and  e).  Kinematic  omega  compares  better 
with  cloudiness;  after  several  hours  of  clearing  in  east  Texas,  clouds 


begin  to  thicken  again  after  1800  GMT  and  thunderstorms  break  out  over 
the  Texas/Louisiana  border  (Fig  6b,  7b,  and  8b). 

F2  omega  ascent  over  much  of  north  and  west  Texas  is  cancelled 
by  FI  omega  descent  (Fig  7d,  f,  and  g). 

Early  Afternoon,  10  April  1979.  A  subsynoptic  surface  low 
(SSL)  develops  around  2100  GMT  (Fig  8a)  and  moves  rapidly  across  the 
Texas  panhandle.  The  clouds  and  weather  (Fig  8b)  change  dramatically: 
a  streak  of  wind-borne  dust  afflicts  west  Texas,  a  thunderstorm  system 
which  spawned  tornadoes  around  Lubbock  now  spreads  over  the  northern 
panhandle  region,  tornadic  thunderstorms  break  out  west  of  Wichita 
Falls,  and  a  comma  cloud  appears  over  Oklahoma.  A  region  of  thun¬ 
derstorms  in  what  has  become  the  “dry  slot"  of  the  comma  cloud,  began 
earlier  southwest  of  Clovis,  New  Mexico,  and  now  showers  are  clearly 
visible  in  satellite  data  above  the  duststorm.  If  one  extrapolates  the 
past  position  of  minor  wave  #2  to  2100  GMT,  it  should  be  in  central 
Kansas  in  Fig  8c.  It  is  very  weak,  and  is  not  illustrated  in  that 
figure.  Minor  wave  #3  is  located  in  north  Texas  above  the  SSL  and  just 
north  of  the  tornado  area  southwest  of  Wichita  Falls,  Texas  (Fig  8a  and 
b).  Minor  wave  #3  is  moving  north-northeastward  at  about  30m/s  over 
the  SSL  moving  east-northeastward  at  9m/s.  Although  positive  vorticity 
advection  would  indicate  pressure  falls  ahead  of  minor  wave  #3,  the 
quasigeostrophic  pressure  tendency  equation  allows  simultaneous  pressure 
rises  there  if  warm  advection  is  stronger  close  to  the  ground  than 
aloft.  The  rapid  development  of  the  SSL  as  minor  wave  #3  passed 
overhead  and  the  outbreak  of  violent  tornadoes  just  in  the  wake  of  the 


minor  wave  suggest  that  dynamical  forcing  associated  with  the  mid-level 
disturbance  Influenced  the  low  level  circulation  very  strongly  for  a 
few  hours.  It  will  be  Interesting  to  see  what  happens  later  as  minor 
wave  #3  gets  farther  away  from  the  SSL  and  warm  front,  and  also  what 
happens  when  minor  waves  #4  and  #5  approach  the  front. 

QG  and  kinematic  omega  (Fig  8d  and  e)  appear  similar. 

However,  QG  omega  shows  near  zero  vertical  motion  on  the  entire  right 
side  of  the  500mb  jet  running  across  western  and  northern  Texas  and 
Oklahoma  (Fig  8c  and  d).  This  compares  poorly  with  the  comma  cloud  and 
thunderstorms  in  Oklahoma  (Fig  8b)  where  kinematic  omega  shows  ascent 
(Fig  8e).  The  reason  for  this  major  failure  of  the  method  is  that  the 
QG  approximation  using  only  FI  and  F2  terms  eliminates  forcing  due  to 
latent  heat  release,  which  is  becoming  important  around  the  convec¬ 
tion  on  the  right  side  of  the  jet,  and  also  ageostrophic  forcing  near 
jet  streaks  (Ucellini  and  Johnson  (1979)).  An  analysis  of  jet  streaks 
and  ageostrophic  accelerations  is  beyond  the  scope  of  this  thesis; 
but  the  reader  is  referred  to  partial  analyses  of  the  same  case  in 
Benjamin  and  Carlson  (1982),  Kocin,  et  al .  (1982),  and  Moore  and 
Fuelberg  (1981). 

Kinematic  omega  (Fig  8g)  indicates  subsidence  in  south  Texas. 
Skies  do  not  clear  in  this  region  (Fig  8b);  in  fact,  one  would  expect 
continued  cloudiness  under  a  strengthening  low-level  inversion. 

However,  a  surface  ridge  is  observed  in  southeast  Texas  beginning  around 
2100  GMT  (Fig  8a)  which  is  associated  with  two  phenomena.  First,  the 
flow  of  moist  air  is  shunted  northwestward  between  Austin  and  San 
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Antonio,  Texas,  so  that  a  tiny  area  northwest  of  Mineral  Wells,  Texas, 
becomes  extremely  unstable  ater  2100  GMT  (isodrosotherm  analysis  not 
Illustrated).  Second,  no  severe  weather  occurred  along  or  south  of  this 
sharp  surface  ridge,  probably  because  of  the  strong  capping  inversion. 

QG  omega  (Fig  8d)  shows  subsidence  just  upwind  of  the 
duststorm  in  west  Texas/New  Mexico  and  rising  motion  with  that  portion 
of  the  comma  cloud  on  the  cyclonic  side  of  the  500mb  jet  (Fig  8c). 
Kinematic  omega  (Fig  8e)  indicates  a  long  area  of  moderate  ascent  over 
the  entire  comma  cloud,  centered  just  east  of  the  most  severe  thun¬ 
derstorms,  and  extending  over  the  "dry  slot,"  including  the  SSL.  In 
summary,  QG  omega  is  poorly  related  to  the  comma  cloud  and  kinematic 
omega  is  well  related  during  the  first  few  hours  of  the  comma  cloud's 
existence. 

FI  and  F2  vary  in  time  and  space  more  than  the  total  QG 
field.  Compare  the  relatively  simple  continuity  of  QG  omega  in  Fig  8d, 
9d,  and  lOd  to  the  interplay  between  FI  and  F2  omega  in  Fig  8f/g,  9f/g, 
and  lOf/g.  The  FI  and  F2  omega  fields  (Fig  8f  and  g)  have  more  maxima 
and  minima  than  the  total  QG  omega  pattern  (Fig  8d).  One  final  obser¬ 
vation  of  interest  is  that  the  FI  omega  pattern  over  Oklahoma  and  north 
Texas  in  Fig  8f  is  upwind  of  the  500mb  vorticity  advection  pattern  of 
Fig  8c,  because  the  trough  axis  Is  vertical  from  the  surface  to  500mb 
above  the  SSL,  but  slopes  southwestward  above  the  mid-troposphere. 

F3  omega  (Fig  8h)  exaggerates  QG  omega  (Fig  8d),  but  also 
shows  a  tiny  subsidence  region  in  north  Texas  just  west  of  the  tornadic 
cells  (Fig  8b). 


savage  tornadoes  have  struck  at  Wichita  Falls,  Texas,  at  Lawton, 
Oklahoma,  and  at  several  towns  in  central  Oklahoma.  The  comma  cloud 
has  tripled  in  area.  The  SSL  is  moving  slowly  and  filling;  the  cold 
front  and  dry  line  move  very  little.  A  line  of  towering  cumuli  deve¬ 
lops  suddenly  east  of  Midland,  Texas,  about  2300  GMT.  Minor  wave  #3  is 
in  northern  Oklahoma  (Fig  9c)  and  minor  wave  #4  appears  in  extreme  west 
Texas.  The  most  severe  thunderstorms  are  occurring  behind  minor  wave 
#3  in  central  Oklahoma  and  just  northeast  of  the  SSL  in  Texas. 

QG  omega  (Fig  9d)  shows  descent  over  the  filling  SSL  behind 
wave  #3.  The  region  of  kinematically  computed  ascent  in  Oklahoma  (Fig 
9e)  is  much  stronger  and  broader  than  the  QG  omega  ascent  area.  The 
most  significant  difference  in  the  two  omega  solutions  is  that  the 
center  of  kinematic  ascent  lies  on  the  500mb  vorticity  ridge  of  minor 
wave  #3  (Fig  9c),  and  extends  south  of  the  wave  slightly  beyond  the  SSL 
and  warm  front  In  Texas.  Figures  7e,  8e,  9e,  and  lOe  Indicate 
excellent  consistency  in  the  kinematic  omega  pattern,  and  the  clouds 
and  weather  at  0000GMT  (Fig  9b)  supports  the  kinematic  analysis. 

Intense  condensation  heating  in  the  northwestern  half  of  Oklahoma  and 
convergence  ahead  of  a  low-level  jet  (not  illustrated)  in  southeast 
Oklahoma/northeast  Texas  are  forcing  vertical  motion  which  QG  omega 
cannot  discern.  In  Fig  9b  and  e,  kinematic  omega  shows  rising  motion 
near  all  active  thunderstorms,  including  the  consistently  observed  "dry 
slot"  thunderstorms  (northwest  of  Amarillo,  Texas),  except  for  the 
upslope  area  around  Colorado  Springs.  Although  the  kinematic  analysis 


shows  ascent  under  the  comma  cloud.  It  shows  moderate  ascent  in  other 
regions  also.  Q6  omega  (Fig  9d)  indicates  ascent  only  under  that  part 
of  the  comma  head  over  and  west  of  the  500mb  jet  (Fig  9b  and  c).  The 
kinematic  method  gives  sharp  definition  to  the  subsidence  pattern  just 
west  of  the  developing  squall  line  near  Midland,  Texas  (Fig  9b  and  e). 

Evening,  10  April  1979.  At  0300  GMT  (Fig  10a  and  b)  the  SSL 
is  nearly  gone.  The  cornua  cloud  doubled  in  size  in  three  hours, 
extending  eastward  as  far  as  Indiana.  Strong  surface  ridging  in  west 
Texas  (Fig  10a)  is  associated  with  development  of  a  squall  line  to  the 
east  (Fig  10b)  which  is  producing  large  hail  and  some  tornadoes. 

Severe  storms  continue  in  central  Oklahoma  east  of  the  SSL  along  the 
warm  front,  but  no  damaging  tornadoes  are  reported  for  a  few  hours. 

Minor  wave  #3  in  southern  Kansas  (Fig  10c)  is  about  300km  past  that 
area  of  weakening  severe  weather.  Minor  wave  #4,  the  cold  front,  and 
the  dry  line  all  intersect  at  the  now  tornadic  squall  line  east  of 
Midland,  Texas. 

The  strongest  QG  omega  ascent  analyzed  in  this  study  (9pnb/s) 
is  in  Kansas  in  a  region  of  little  or  no  convective  activity  (Fig  10b 
and  d).  It  encompasses  the  part  of  the  comma  head  west  of  the  500mb 
jet  (Fig  10c),  but  it  also  invades  the  dry  slot  where  one  might  expect 
subsidence.  Kinematic  omega  (Fig  lOe)  shows  most  of  the  comma  in  a 
region  of  500mb  ascent;  however,  it  also  shows  moderate  dry  slot 
ascent,  similar  to  the  QG  solution.  The  coma  extends  northeastward 
over  a  region  of  subsidence  In  Iowa.  The  QG  solution  does  not  show  much 
ascent  over  the  strong  convection  in  southwest  Missouri,  although  kine- 


matlc  omega  does.  Both  solutions  show  moderate  or  better  ascent 
associated  with  the  minor  wave  In  southern  Kansas  with  an  Important 
difference:  QG  omega  shows  ascent  changing  to  descent  as  the  upper 
wave  passes,  while  kinematic  omega  shows  strong  ascent  centered  on 
the  troughline  of  the  minor  wave,  extending  several  hundred  kilometers 
upstream  of  the  trough.  QG  omega  shows  subsidence  behind  wave  #4  In 
west  Texas,  but  kinematic  omega  shows  rising  motion  there  (Fig  10c, 
d,  and  e).  The  surface  ridge  In  south  Texas  (Fig  10a)  shifts  south¬ 
ward  under  continued  (kinematically  computed)  subsidence. 

Beginning  at  this  time,  QG  omega  (Fig  lOd)  shows  Intense 
subsidence  over  west  Texas  and  southern  New  Mexico  In  contradiction 
to  kinematic  omega  (Fig  lOe).  Dlabatlc  cooling  of  the  surface  would 
Intensify  subsidence,  so  some  other  ageostrophlc  forcing  must  account 
for  this  unreasonable  QG  result,  the  area  Is  also  the  scene  of  moder¬ 
ate  F4  (deformation-induced)  forcing.  In  Krlshnamurtl 's  (1968)  case 
study,  two  terms  compensated  for  very  strong  forcing  due  to  thermal 
advectlon  In  a  synoptic  location  very  similar  to  his  case  (see  his 
Figures  10  and  11).  These  were  the  "differential  deformation"  and 
"differential  divergence"  terms  In  the  balance  omega  equation.  The 
location  was  just  upstream  from  an  occluding  frontal  wave,  equatorward 
from  the  500mb  low  (and  slightly  upstream),  on  the  left  flank  of  the 
500mb  jet.  These  two  ageostrophlc  effects  may  have  compensated  for 
strong  cold-advectlon  subsidence  on  11  April  1979  as  well. 

F3  omega  (Fig  lOf)  begins  to  have  significant  deviations 
from  QG  omega  (Fig  lOh)  at  this  time.  In  Its  analysis  In  west  Texas, 
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Kansas,  and  Oklahoma,  QG  omega  complements  the  cloud  and  weather  pat¬ 
terns  (Fig  10b)  less  and  less.  So  this  failure  of  Wi in-Nielsen's 
method  to  approximate  QG  omega  comes  at  a  time  when  QG  omega  itself 
looks  rather  poor. 

Midnight,  10-11  April  1979.  The  SSL  is  gone  by  0600  GMT  (Fig 
11a  and  b),  and  the  cold  front  slows  down  in  central  Texas.  The  squall 
line  in  north  Texas  becomes  coincident  with  the  tail  of  the  comma  cloud. 
Thunderstorms  persist  in  the  dry  slot  of  the  comma  cloud.  Three  minor 
waves  are  evident  in  the  500mb  relative  vorticity  pattern  (Fig  11c). 

Wave  #3  is  on  the  Kansas/Nebraska  border,  wave  #4  is  in  western 
Oklahoma,  and  wave  #5  is  in  southeastern  New  Mexico.  Although  some 
severe  weather  continues  in  northeastern  Oklahoma,  the  most  severe  thun¬ 
derstorms  are  occurring  at  the  Intersection  of  minor  wave  #4  and  the 
warm  front, and  behind  wave  #4  along  the  cold  front  in  central  Texas. 

The  QG  and  kinematic  omega  fields  are  not  even  remotely  simi¬ 
lar  at  this  time  (Fig  lid  and  e).  The  QG  pattern  is  generally  perpen¬ 
dicular  to  the  comma  cloud  tail,  and  some  of  the  active  thunderstorms 
are  in  areas  of  QG  omega  subsidence  (Fig  lib).  The  kinematic  omega 
pattern  is  very  similar  to  the  comma  cloud  except  that  (1)  the 
strongest  kinematic  ascent  is  along  the  poleward  side  of  the  comma  tail. 


and  (2)  the  part  of  the  comma  head  to  the  left  of  the  500mb  jet  (Fig 
11c)  is  indicated  as  subsiding.  Thus,  the  center  line  of  the  comma 
tail  is  slightly  downwind  (northeast)  of  the  axis  of  maximum 
kinematically-computed  ascent.  The  strongest  ascent  of  the  entire 
study  (13*kb/s)  is  analyzed  200km  northeast  of  the  nose  of  the  dry  slot. 


Strong  Q6  omega  descent  (Fig  lid)  In  New  Mexico  is  largely 
the  result  of  F2  omega  (Fig  Ilf),  but  FI  omega  (Fig  llg)  shows  strong 
rising  motion  in  north  Texas. 

Early  morning,  11  April  1979,  A  mesolow  forms  on  the  cold 
front  near  the  Texas/Oklahoma  border  at  0900  GMT  (Fig  12a).  Both  minor 
wave  #4  in  northern  Kansas  and  minor  Wave  #5  in  west  Texas  are  too 
distant  from  this  region  to  account  for  this  development  (Fig  12c). 

The  squall  line  in  central  Texas  is  moving  northeastward,  well  north  of 
the  surface  ridge  on  the  Texas  coastal  plains  (Fig  12a  and  b).  Less 
severe  weather  in  eastern  Oklahoma  and  north  Texas  was  reported  at  this 
time  of  the  morning  than  in  earlier  hours. 

Kinematic  omega  (Fig  12e)  conforms  to  the  comma  cloud  pattern 
very  well  (Fig  12b),  except  that  the  ascent  In  Texas  is  much  more 
widespread  than  the  clouds  and  thunderstorms.  There  has  been  a  con¬ 
sistent  pattern  since  the  comma  cloud  first  appeared  around  2100  GMT 
(see  Fig  8b/e,  9b/e,  lOb/e,  lib/e,  and  13b/e)  for  the  kinematic  omega 
to  Indicate  ascent  south  and  southeast  of  the  tip  of  the  comma 
cloud,  and  for  the  comma  cloud  tail  to  grow  toward  the  southwest  with 
the  repeated  squall  lines  in  north  Texas.  Most  of  the  comma  cloud 
advects  downstream  very  rapidly,  but  the  tip  of  the  comma  cloud  tail 
stays  in  north  Texas  in  a  region  of  rising  motion.  QG  omega  (Fig  12d) 
shows  descent  just  northeast  of  an  active  squall  line  in  north  Texas 
and  ascent  over  all  of  the  dry  slot  (Fig  12b).  There  j_s  a  thunderstorm 
in  the  dry  slot,  but  the  kinematically  computed  (Fig  12e)  descent  over 
most  of  the  dry  slot  (except  zero  at  the  thunderstorm  location)  looks 
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more  believable.  In  northern  New  Mexico,  downslope  winds  (see  the  sur¬ 
face  omega  In  Fig  12a)  are  reflected  in  the  kinematic  omega  but  not  in 
QG  omega. 

FI  and  F2  omegas  (Fig  12f  and  g)  both  contribute  to  extreme 
QG  omega  subsidence  in  southeastern  New  Mexico  (Fig  12d).  F3  omega 
(Fig  12h)  seems  to  moderate  the  QG  omega  maximum.  Recall  that  F3  omega 
normally  overestimates  QG  maxima.  F4  omega  (Fig  12i)  rises  to  +4^.b/s 
near  a  SSxlO1^-!  relative  vorticity  maximum  in  New  Mexico. 

Dawn,  11  April  1979.  By  1200  GMT  (Fig  13a  and  b)  most  of  the 
comma  cloud  has  moved  out  of  the  analysis  region.  Very  little  severe 
convection  continues.  However,  minor  wave  #5  (Fig  13c)  is  in  central 
Texas  approaching  the  mesolow  on  the  cold  front.  At  this  point,  one 
might  forecast  renewed  severe  weather  activity  east  and  northeast  of 
this  mesolow  beginning  about  the  time  wave  #5  would  be  above  the  meso¬ 
low  (say,  1830  GMT).  Severe  weather  did  break  out  again  at  1850  GMT 
south-southeast  of  the  mesolow,  and  it  spread  northeastward  into 
Oklahoma,  Louisiana,  and  Arkansas  during  11  April  1979  (Alberty,  £t 
iU  (1979) ). 

Again,  QG  omega  (Fig  13d)  seems  unrelated  to  the  comma 
cloud,  thunderstorm,  or  kinematic  omega  patterns.  Kinematic  omega 
(Fig  13e)  compares  very  favorably  with  clouds  and  weather  (Fig  13b). 

F3  omega  (Fig  13h)  exaggerates  QG  omega  (Fig  13d)  in  an  ana¬ 
lyzed  ascent  area  in  central  New  Mexico.  A  maximum  of  relative  vor¬ 
ticity  (Fig  13c)  dominates  the  F4  omega  pattern  (Fig  13i). 


B.  Richardson  Humber 


In  areas  of  low  Richardson  number  (Rj),  turbulent  motion  Is 
favored  more  than  In  other  areas.  This  motion  moves  heat,  moisture,  and 
momentum  upward  from  the  surface  much  more  efficiently  than  the  rate  by 
molecular  diffusion.  If  areas  of  low  R-j  exist  near  the  jet  stream  level, 
strong  horizontal  momentum  can  be  exchanged  downward.  For  these  reasons, 
R-j  was  analyzed  independently  in  two  vertical  regimes,  one  below  and 
one  above  500mb.  Two  criteria  for  low  values  of  Rj  are  suggested  by 
Stone  (1966).  0.25<R1-<.95  may  be  associated  with  symmetric  instability 
and  0<R.j < .25  may  be  associated  with  Kelvin-Helmholz  instability. 

Figures  14a/b  through  22a/b  portray  general  areas  of  low 
Richardson  number  in  two  vertical  regimes:  mid/upper  troposphere  ( 500- 
lOOmb)  and  boundary  layer/lower- troposphere  (near  surface  to  500mb). 

Each  figure  shows  the  region  with  0<R^<.95  and,  enclosed  within  this 
area,  with  0<R^<.25.  Considering  Rj  as  a  continuous  mathematical 
function  of  height,  it  is  impossible  to  identify  a  sounding  as  having 
some  layer  with  0<Rf<.25  unless  the  neighboring  layer  above  and  below 
has  .25<R-j<.95.  That  is,  a  sounding  with  a  layer  with  the  smaller 
values  of  Rj  must  have,  therefore,  some  nearby  layer  with  the  larger 
values  (0<R^< .95)  as  well.  Areas  with  . 25<R^<.95  are  easy  to 
identify  in  Fig  14a/b  through  22a/b  between  the  outer  (solid)  contour 
and  the  inner  (dashed)  contour.  In  Addition,  Fig  14b,  15b,  16b,  etc. 
illustrate  areas  where  R^.95  in  the  lowest  25mb  above  the  surface, 
including  reports  of  negative  lapse  rates  of  potential  temperature. 
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The  fields  of  Richardson  number  suggest  widespread 
instability  with  details  that  vary  considerably  over  three-hour 
periods.  Fig  14a  and  15a  compared  show  that  low  values  of  Ri  above 
500mb  spread  behind  minor  wave  #1  crossing  Kansas.  From  1800  to  2100 
GMT  (Fig  16a  and  17a)  minor  wave  #3  leaves  a  region  of  low  R^  in  its 
wake  across  west  Texas,  and  has  a  relatively  stable  region  in  central 
Texas  to  its  south.  This  instability  continues  following  minor  wave  #3 
through  0300  GMT  (Fig  18a  and  19a)  across  Oklahoma.  Minor  wave  #4 
exhibits  a  very  similar  pattern  as  it  crosses  northwest  Texas  and 
western  Oklahoma  in  Fig  19a  through  22a.  Low,  even  very  low  Richardson 
numbers  spread  on  the  poleward  side  of  its  wake,  and  a  curious  stable 
spot  occurs  south  or  southeast  of  the  minor  wave  trough.  By  1200  GMT 
(Fig  22a)  two  of  these  stable  regions  (one  in  Arkansas,  one  in  central 
Texas)  and  two  regions  of  instability  (one  in  Kansas,  one  in  west 
Texas)  seem  to  be  associated  with  minor  waves  #4  (in  Nebraska)  and  #5 
(in  southwest  Oklahoma). 

Below  500mb  (Fig  14b-22b),  low  Richardson  numbers  are  found 
not  only  in  the  area  behind  the  dry  line,  but  also  around  and  often 
northeast  of  the  mesolow  in  north  Texas,  and  in  the  warm  sector 
southeast  of  the  mesolow.  The  unstable  region  around  the  mesolow 
increases  in  area  three  times  (Fig  17b,  19b,  and  22b),  possibly  in 
phase  with  the  500mb  minor  waves  (compare  to  Fig  17a,  19a,  and  22a).  A 
possible  explanation  for  this  phenomenon  might  be  that  transient  wind 
maxima  associated  with  the  minor  waves  cause  ageostrophic  low-level 
accelerations  northwestward  under  the  exit  region  of  the  upper  tro- 


pospherlc  jet.  With  each  pulse  of  the  low  level,  warm  flow  under  the 
upper  jet  (see  the  85Qmb  analyses  in  Moore  and  Fuelberg  (1981)),  "b&ftz 
decreases  and  3|  u  )  f)z  increases  momentarily. 

Figure  16b  shows  such  a  region  around  the  Red  River  at  the 
time  severe  storms  are  breaking  out  in  north  Texas,  even  before  the  SSL 
arrives. 


CHAPTER  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions.  Five  mid-level  disturbances  of  wavelength  500  to 
800  km  passed  through  the  area  at  speeds  of  20  to  33  m/s,  evident  in 
strong  variations  of  mid-level  height,  vorticity,  vertical  motion,  and 
Richardson  number  on  a  three-hour  scale.  The  500mb  quasi-geostrophic 
(QG)  omega  field  often  showed  rising  followed  by  subsidence  as  these 
disturbances  passed,  while  kinematically  computed  omega  showed 
ascent  on  the  troughline  of  these  disturbances.  The  discrepancy  between 
the  two  omega  solutions  was  probably  a  result  of  physical  forcing  not 
considered  in  the  simplified  QG  omega  equation,  namely  latent  heat 
release  and  ageostrophic  response  in  general.  The  "minor"  short  waves 
were  associated  with  movement  of  the  subsynoptic  surface  low  (SSL)  and 
dry  line,  with  development  and  filling  of  the  SSL,  with  the  location  and 
onset  of  severe  storms,  and  with  major  features  of  the  comma  cloud.  As 
two  of  these  waves  passed  the  SSL  or  the  cold  front/dry  line  intersec¬ 
tion  in  central  Texas,  severe  thunderstorms  broke  out  or  intensified 
just  northeast  through  southeast  of  the  surface  low  and  remained  violent 
until  the  waves  were  several  hundred  kilometers  downstream. 
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A  comparison  of  vertical  motion  computed  by  the  kinematic 
method  (with  O'Brien  correction)  and  the  QG  omega  equation  method,  to 
clouds  and  weather  indicates  that  both  methods  are  able  to  resolve 
detail  in  weather  phenomena  as  small  as  three  gridlengths  (480km)  under 
the  best  circumstances.  However,  QG  omega  became  very  inaccurate  as 
ageostrophic  processes  became  stronger  during  this  study.  In  fact,  both 
methods  give  spurious  results  at  times  and  one  can  make  the  following 
statement:  for  research  purposes  these  methods  make  an  independent  and 
complementary  pair,  while  for  operational  purposes  the  computational 
ease  and  ageostrophic  sensitivity  of  the  kinematic  method  make  it 
preferable  to  the  QG  method. 

In  comparison  of  the  QG  omega  and  kinematic  omega  to  clouds 
and  weather,  the  QG  method  apparently  failed  in  many  situations.  It 
grossly  overestimated  the  subsidence  region  in  southern  New  Mexico 
southwest  of  the  closing  upper  low.  Krishnamurti 's  (1968)  solution  of 
the  balance  wind  omega  equation  for  a  very  similar  synoptic  situation 
showed  an  almost  identical  subsidence  region  forced  by  cold  advection, 
which  was  counterbalanced  by  two  ageostrophic  forcing  functions:  dif¬ 
ferential  deformation  and  differential  divergence.  Kinematic  omega  con¬ 
sistently  analyzed  rising  motion  in  active  thunderstorm  areas.  Its  abi¬ 
lity  to  locate  ascent  forced  by  physical  factors  neglected  in  the  QG 
omega  equation  which  was  solved,  especially  latent  heat  release  and 
ageostrophic  accelerations  under  the  exit  region  of  a  jet  streak, 
establish  the  kinematic  method  as  definitely  superior  to  the  QG  method 
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for  forecasting  purposes.  Sometimes  kinematic  omega  indicated  ascent 
In  an  area  before  thunderstorms  occurred. 


The  QG  method  consistently  analyzed  500mb  ascent  under  that 
portion  of  the  comma  "head"  to  the  left  of  the  jet,  but  elsewhere 
failed.  In  general  ,  the  QG  analysis  did  not  recognize  significant 
vertical  motion  on  the  anticyclonic  side  of  the  500mb  jet.  One  may 
infer,  tentatively,  that  vertical  motion  due  to  latent  heat  release 
and  to  ageostrophic  forcing  on  the  right  of  the  mid-tropospheric  jet 
became  as  significant  as  the  quasi -geostrophic  forces.  In  particular, 
other  researchers  have  found  evidence  that  the  accelerations  in  the 
exit  region  of  a  jet  streak  became  important  (Kocin,  et  al .  (1982), 
Moore  and  Fuelberg  (1981)).  As  shown  in  Ucellini  and  Johnson  (1979), 
strong  vertical  motion  occurs  in  the  low  levels  of  jet  streak  exit 
regions,  induced  by  ageostrophic  forces. 

The  relationship  of  kinematically  computed  omega  to  the 
comma  cloud  can  be  summarized  as  follows. 

1.  Kinematic  omega  was  always  centered  on  the  active  convec 
tion  on  the  western  edge  of  the  comma's  tail,  so  that  the  center  line 
of  the  comma  tail  advected  in  time  slightly  east  of  the  axis  of 
maximum  instantaneous  rising  motion. 

2.  New  centers  of  ascent  formed  repeatedly  southeast  and 
south  of  the  comma  tall  before  the  tall  extended  over  these  new 
convective  regions. 

3.  Kinematic  calculations  indicated  that  omega  In  the  "dry 
slot"  varied  considerably,  but  tended  to  be  upward  around  the  con¬ 
sistently  observed  dry  slot  thunderstorms.  The  sharp  back  edge  of  the 


comma  cloud  can  not  be  taken  as  a  dividing  line  between  Instantaneous 
descent  and  ascent. 

4.  Kinematic  omega  usually  Indicated  rising  motion  under 
that  portion  of  the  comma  head  to  the  left  of  the  500mb  jet;  however. 

In  one  analysis  (0600  GMT)  the  opposite  was  true. 

Two  details  of  this  analysis  apply  specifically  to  forecasting 
severe  thunderstorms  In  the  southern  plains  In  springtime.  First,  both 
the  kinematic  and  the  QG  analyses  showed  an  area  of  500mb  rising  motion 
in  eastern  New  Mexico  at  least  six  hours  before  an  SSL  and  tornadoes 
formed  just  to  the  east  In  Texas.  Second,  subsidence  over  the  Texas 
coastal  plain  caused  ridging  and  the  maintenance  of  an  Intense  low- 
level  inversion.  These  features  simultaneously  surpressed  convection 
In  southeast  Texas  and  diverted  moisture  northwest  of  the  region. 

Minor  waves  aloft  brought  on  very  low  Richardson  number  (Rj) 

In  the  middle>  and  high-troposphere,  plus  small  regions  of  high  R ^  to 
their  south.  Below  500mb,  very  low  Rf  was  typical  of  the  area  west  of 
the  dry  line,  the  SSL,  the  warm  sector  southeast  of  the  SSI,  and  some¬ 
times  the  warm  front  area  northeast  of  the  SSL.  Rj  analyses  showed 
strong  vertical  wind  speed  shear  directly  over  the  SSL. 

Recommendations  for  further  study.  Many  mechanisms  In  this 
storm  are  contained  In  Its  three-dimensional  structures  of  horizontal 
wind,  vertical  motion,  and  other  parameters.  A  complete  analysis, 
particularly  a  three-dimensional  analysis  of  moisture,  clouds,  and 
weather,  could  provide  a  definitive  portrait  of  the  evolution  of  the 
system.  This  thesis  highlights  three  points  of  departure.  First,  more 

Information  about  the  structure  of  meso-a  scale  upper  disturbances 
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which  are  termed  "minor  waves"  in  this  thesis  and  their  relationship  to 
the  low  level  circulation  and  weather  is  necessary  to  the  understanding 
of  the  organization  of  cloud  patterns  and  severe  weather.  The  problem 
of  detecting  and  forecasting  these  waves  is  important.  Second,  the 
implications  of  ageostrophic  accelerations  In  the  exit  region  of  jet 
streaks  need  to  be  explored  and  modeled  on  smaller  scales  (meso-a  or 
meso-B)  than  studies  heretofore  (Ucellini  and  Johnson  (1979)). 

Third,  a  fully  integrated  moisture  and  cloud  analysis  every  three  hours 
would  be  an  ideal  verification  tool  for  all  types  of  cloud  forecasting 
techniques. 

A  study  of  the  magnitude  and  source  of  the  subsidence  in 
south  Texas  may  be  able  to  quantify  the  contributions  of  adiabatic 
descent  and  thermal  advection  toward  maintenance  of  the  capping  inver¬ 
sion  In  that  region.  Carlson  ,  et  al .  (1980)  have  proposed  that  the 
advection  of  warm  air  heated  over  the  Mexican  Plateau  Is  the  dominant 
process  that  maintains  this  Inversion. 

Calculations  were  made  In  a  pressure  coordinate  system  rather 
than  a  o-coordinate  system  in  hopes  of  achieving  more  consistent 
vertical  finite  differencing.  Unfortunately,  (1)  a  number  of  grid- 
points  below  ground  level  were  rendered  useless  in  the  process, 

(2)  the  computer  run  time  required  to  test  for  flagged  grldpoints  is 
excessive,  and  (3)  the  task  of  programming  a  computer  to  test  for 
flags  on  each  use  of  grldded  data  proved  to  be  very  complicated  and 
time-consuming.  Oespite  the  straightforward  simplicity  and  apparently 
more  consistent  vertical  differencing,  the  pressure  coordinate  system 
is  less  practical  than  the  o-coordinate  system. 
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Preparation  of  Data  Before  Interpolation  to  a  Grid 


The  AVE-SESAME  1  case  contains  a  total  of  335  rawinsonde 
ascents  made  from  39  stations  every  three  hours  from  1200  GMT  10  April 
1979  through  1200  GMT  11  April  1979.  Figure  1  illustrates  the  geo¬ 
graphic  location  of  these  stations.  Data  were  collected  at  all  stations 
and  at  all  times  except  those  indicated  as  missing  in  Table  6. 


Missing 

"Questionable 

1200 

10  April 

CLL,  FSM 

BVD, 

GAG,  GLD 

1500 

10  April 

CLL 

GAG 

1800 

10  April 

CLL 

BVO, 

GAG 

2100 

10  April 

MRF 

FSM, 

GAG 

0000 

11  April 

BVO,  MTX 

SEP, 

GAG 

0300 

11  April 

None 

GAG 

0600 

11  April 

ABI,  CLL, 

MTX 

LCH, 

GAG 

0900 

11  April 

BVO,  COU, 

RTN 

GAG 

1200 

11  April 

BVO,  DUA, 

UOX 

FSM, 

GAG 

Table  6:  Missing  rawinsondes  and  rawinsondes  which  the  NASA  archivers 
considered  questionable. 


Although  NASA  archivers  considered  sixteen  of  the  rawinsondes 
"questionable",  these  data  were  used  without  reservation  in  this  analysis. 
A  subjective  analysis  explained  in  Section  B  of  Chapter  II  indicated  the 
following  biases  in  height  data  from  four  stations:  GAG(+17m),  ABI(+12m), 
SEP(+7m),  and  GLD(-5m).  All  height  reports  from  these  stations  were 
corrected  before  the  objective  analysis. 


Alt  data  at  all  pressure  levels  were  interpolated  linearly  to 
a  balloon  release  one  hour  before  data  time.  For  example,  information 
from  rawinsondes  released  at  Abilene,  Texas,  at  1121,  1442,  and  1740  GMT 
were  used  to  construct  a  sounding  which  would  have  resulted  from  a  balloon 
released  at  1400  GMT,  and  these  data  bear  the  label  "1500  GMT  sounding". 
More  explicitly,  the  linear  interpolation  to  1400  between  1121  and  1442 
was  averaged  with  the  linear  extrapolation  from  1740  past  1442  to  1400. 

Surface  pressure  tendencies  were  calculated  at  each  time  at 
each  station  from  the  rawinsonde  reported  surface  pressure.  Because  no 
data  was  available  before  the  first  data  set  at  1200  GMT,  10  April  1979, 
the  computed  surface  pressure  change  designated  for  that  time  is  actually 
valid  at  1330  GMT.  Likewise,  the  1200  GMT,  11  April  1979  change  is 
valid  at  1030  GMT. 

Some  data  for  the  12-hourly  standard  rawinsonde  were  computed  by 
hand.  Height,  wind  direction,  and  windspeed  at  non-mandatory  levels  (1000, 
900,  800,  and  600mb)  were  calculated  by  logarithmic  interpolation  of 
Service  C  data(p3<P2<Pj): 

*2*zl  *  (Rd/g)  Tv(l,2)ln(Pl/Pz)  (A.l) 

Z3-Z1  =  (Rd/g)  Ty(i,3)ln(pi/P3)  (A. 2) 

where  P2  refers  to  a  pressure  level  between  two  standard  pressure  sur- 
4  faces  pj  and  P3,  while  Zi  and  Z3  are  reported  heights.  Rd  is  the  gas 

constant  and  g  is  the  acceleration  of  gravity.  Notice  that  the  highest 
pressure,  pj,  may  be  the  station  pressure  if  Z 3  is  taken  as  the  station 
4  elevation.  Now,  If  the  mean  virtual  temperature  Tyj 1 >2)  in  the  P1-P2  layer 

is  approximately  equal  to  Tv( 1^3)  in  the  P1-P3  layer,  the  equations  can  be 
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(A. 3) 


combined  to  yield  the  height  Z 2  at  P2 : 

Z2=Z1+(Z3-Z1)ln(p1/p2)/ln(p1/p3) 

The  error  introduced  by  assuming  Tv  constant  causes  Z2  to  be  estimated 
too  small,  because  Tvq  3)  is  less  thatTv(ij2)  generally  (recall  P3<P2<Pl ) 
One  may  estimate  the  error  at  6Q0mb,  where  it  is  largest  in  absolute 
value  due  to  the  logarithmic  spreading  of  pressure  layers  with  height. 
Assuming  TV=T,  the  solution  of  Equation  (A. 3)  for  NACA  Standard 
Atmospheric  heights  (Haltiner  and  Martin  (1957))  for  P3=500mb,  p2=600mb, 
Pl*700mb  yields  a  600mb  height  estimate  of  4183m.  This  is  17m  less 
than  the  standard  atmosphere  4200m.  This  research  did  not  make  an  adjust¬ 
ment  for  the  expected  error  in  the  height  estimate,  but  simply  used  the 
result  of  Equation  (A. 3). 

Terrain  heights  were  derived  from  two  data  sources  and  then 
smoothed.  For  gridpoints  on  or  east  of  105°W,  a  U.S.  Air  Force  data  set 
(Carr  (1981))  at  1°  latitude-longitude  Interval  was  interpolated  linearly, 
while  west  of  105°W,  Heights  were  estimated  from  a  1:3  million  scale 
weather  chart. 


Appendix  B 

Gridded  Analysis  Calculations 

The  grid  chosen  for  this  study  is  an  11x12  array  on  a  polar 
stereographic  projection,  with  grid  spacing  in  both  directions  one  inch 
at  1:7.5  million  scale.  The  reference  meridian  is  100°  west  longitude 
(Figure  1).  At  35°  north  latitude,  161km  separates  grid  points,  and 
this  is  taken  as  the  nominal  grid  distance.  The  average  distance  of 
each  station  from  its  nearest  neighbor  for  SESAME  stations  inside  the 
grid  (i.e.  west  of  91°  west)  is  214km,  or  1.33  gridpoints.  The  grid 
extends  vertically  every  lOOmb  from  lOOOmb  to  lOOmb,  plus  a  surface  grid. 

Cressman  analysis.  Data  were  interpolated  to  the  grid  using  a 
Cressman  (1959)  weighting  function: 

R2  -  d2 

w  =  — - -  for  <kR  (B.l) 

r2  ♦  d2 

where  R  is  a  constant  radius  of  influence  and  d  is  the  distance 
between  the  gridpolnt  and  the  balloon  location  at  the  given  pressure 
level.  Figure  23  shows  that  the  weighting  function  falls  from  one 
at  the  grid  point  to  zero  at  and  beyond  distance  R. 
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Only  data  within  the  radius  of  influence  are  used  to  calculate 


grid  values  according  to  the  formula: 

n  n 

S0  *  (EwiSi)/Zwi  ( B . 2 ) 

i=l  1=1 

where  data,  ,  observed  at  locations  (x^.y^)  within  the  radius  of 
influence,  R,  of  a  particular  gridpoint  (x0,  y0)  are  weighted  relative 
to  the  distances  from  the  gridpoint. 

The  height  analysis  utilizes  the  observed  wind  at  each  data  point 
to  estimate  the  height  gradient.  After  eliminating  second  order  terms 
of  the  Taylor's  series  expansion  of  Z(x,y), 

3Z  dZ 

Z(x+ax,y+ay)  *  Z{x,y)  +  —  ax+  —  a y  (B.3) 

3  x  3y 

But  the  observed  wind  should  be  a  good  approximation  of  the  geostrophic 
wind  (Cressman,  1959),  except  that  the  observed  wind  is  typically 
weaker  than  the  geostrophic  wind  by  a  factor  of  1.08: 

3Z 

—  =  (f0/g)v0bs  x  1-08 

3X 

3Z 

—  *  (-^o/9^uobs  x  1<08  (B.4) 

3y 

for  f0  *  constant  value  coriolis  parameter  and  g  =  acceleration  of  gravity. 

The  equation  used  to  calculate  height  at  gridpoints  was 

n  3Z  3Z  n 

wiki ( Z^+ — (xf-x0)+— (yi-y0})/E(wfkf)  (B.5) 

1*1  3x  3y  1*1 


with  kj*l  when  no  wind  data  were  available  and  kj=4  otherwise. 
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After  the  second  pass  approximation  to  the  grldded  fields,  miss¬ 
ing  height  and  Mind  data  were  bogused  from  the  latest  approximation,  and 
these  bogus  data  were  treated  equally  with  observed  data  In  the  next 
analysis  pass.  Thus,  there  were  no  missing  data  In  the  working  rawin- 
sonde  array  after  the  second  pass. 

Missing  data  sometimes  made  it  necessary  to  increase  the  fixed 
radius  of  influence  in  increments  of  0.5  gridlengths  at  individual  points 
until  at  least  two  data  points  were  included.  But  a  listing  of  radii 
greater  than  3  gridlengths  showed  that  R=5.0  was  the  maximum  radius  (on 
the  extreme  northwest  and  southwest  corners  of  the  grid).  Moreover, 
when  12-hourly  data  were  available,  R  was  nowhere  larger  than  3.5. 

Smoothing.  The  simple  smoothing  technique  mentioned  throughout 
this  thesis  involves  the  double  application  of  two  one-dimensional 
operator s: 

1  s 

T1,j  =  <  l-s)fi ,  j  +  2(f1+l,j+fi-l,j)  (B.6) 

j  s 

T1,j  =  d'S^iJ  ♦  2(fi,j+l+f1,j-l>  (B-7) 

First  (B.6)  is  applied  to  the  entire  field,  with  s=+£,  then  (B.7)  is 
applied  with  s=+^.  Next  (B.6)  is  applied  with  s=~;  and  finally  (B.7) 
with  s*~.  The  response  function 

R(lx,ly)  -  ( 1-si n4rr d/1  x)(l-s1n4nd/ly)  (B.8) 

shows  the  amplitude  of  smoothed  disturbances  of  wavelength  lx  and  ly 
in  the  x  and  y  directions. 

T(lXfly)  *  f  x  R(lx,1y)  (B.9) 

In  particular  2d  waves  (in  either  direction)  are  completely  quashed. 
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Bilinear  Interpolation.  Estimated  values  of  height  and  wind 
components  were  derived  by  bilinear  interpolation  of  gridded  data  to  the 
location  of  the  rawinsonde.  With  reference  to  Figure  24  ,  the  value  of 
field  S  at  point  P5  between  gridpoints  g^  and  g2  is  given  by 

S5  =  S^d-a)  +  S2x(a)  (B.10) 

and  at  point  P6  by 

Sg  *  S3x(l-a)  +  S4XU)  (B.ll) 

Then  the  value  of  S  at  the  target  p0  between  P5  and  pg  is  given  by 
S0  =  S5x( 1-b)  +  S6(b)  (B.12) 

This  same  method  is  used  later  to  interpolate  in  the  vertical  and  time 
domains  for  the  trajectory  calculations  (trl linear  and  quadri linear 
interpolation) . 

Error  calculations.  Root-mean-square  (r.m.s.)  deviations  of 
the  SESAME  data  set  were  calculated  by 
2  9  Nt  11  12  2 

(oz(p)J  =  l  ll  (I  l  (z0(P>-Za(P»l  J)/NZ  ( B .  13 ) 

t=l  n=l  i=l  j=l  tP 

and 

o  9  Nt  11  12  . 

(°v(p)  )2  *  l  (I  (I  l  C(u0(p)-ua(p))2  (B.14) 

t=l  n=l  i=l  j=l 

+  (v0(p)-va(p))Z]tp)]/Nv 

for  each  pressure  level.  refers  to  the  number  of  soundings 
at  time  t  (1=12Z/10  April,  2=15Z/10  April,  etc.).  The  notation  |tp 
indicates  that  each  piece  of  sounding  data  was  compared  to  the  analysis 
for  the  same  time  and  pressure  level.  Nz  and  Nv  refer  to  the  total 
number  of  deviations  that  were  computed,  discounting  missing  data. 
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time (GMT) 

1200 

1500 

1800 

2100 

0000 

0300 

0600 

0900 

1200 

°z 

pressure(mb) 

100 

20.2 

15.0 

11.4 

7.7 

11.5 

33.0 

17.2 

17.3 

14.4 

24.7 

200 

9.6 

9.8 

12.1 

39.8 

23.8 

23.1 

16.6 

14.0 

10.5 

30.5 

300 

8.8 

9.9 

10.7 

23.3 

27.4 

13.6 

17.6 

12.1 

12.1 

28.5 

400 

6.5 

7.7 

7.4 

9.7 

11.4 

7.6 

10.3 

7.9 

9.3 

22.6 

500 

5.5 

6.6 

5.5 

6.5 

8.8 

6.9 

6.4 

10.0 

11.4 

19.6 

600 

6.8 

4.6 

4.7 

6.2 

5.9 

6.8 

6.7 

7.8 

7.2 

18.8 

700 

6.2 

4.1 

5.0 

7.1 

4.6 

5.3 

5.6 

6.0 

8.0 

15.1 

800 

7.9 

3.8 

5.2 

5.2 

4.9 

5.2 

5.2 

7.2 

5.5 

16.9 

900 

11.2 

4.2 

4.7 

4.9 

5.0 

7.9 

11.3 

5.8 

4.4 

14.8 

1000 

1.5 

2.3 

2.2 

.9 

1.4 

2.5 

2.6 

.2 

1.4 

8.5 

Table  7: 

R.m.s. 

height  error  (m) 

with  R; 

=2d  after  five  passes. 

time (GMT) 

1200 

1500 

1800 

2100 

0000 

0300 

0600 

0900 

1200 

°v 

pressure(mb) 

100 

1.4 

4.5 

1.6 

2.6 

2.1 

2.9 

9.4 

13.8 

7.7 

6.5 

200 

5.3 

6.3 

4.5 

5.3 

3.8 

4.8 

9.1 

13.2 

6.5 

9.2 

300 

3.9 

5.7 

4.3 

8.5 

5.8 

6.1 

8.8 

6.8 

7.7 

11.3 

400 

3.4 

3.7 

3.8 

6.2 

6.6 

6.0 

5.9 

4.8 

4.7 

9.8 

500 

3.1 

3.2 

3.4 

5.0 

4.9 

4.8 

8.2 

6.4 

4.0 

8.8 

600 

2.6 

3.7 

3.0 

3.1 

5.4 

4.4 

7.4 

4.7 

3.2 

7.3 

700 

5.0 

2.4 

2.4 

2.8 

6.3 

4.6 

5.9 

4.8 

5.1 

6.5 

800 

2.5 

2.5 

2.5 

2.5 

5.3 

5.4 

5.6 

4.1 

4.6 

6.6 

900 

2.0 

2.2 

2.1 

2.1 

3.0 

3.6 

4.6 

4.5 

3.3 

5.4 

1000 

.9 

1.2 

.7 

.7 

.8 

2.4 

2.0 

.8 

2.0 

3.6 

surface 

2.1 

2.6 

2.5 

3.0 

3.3 

4.2 

3.8 

2.9 

2.1 

4.8 

Table  8:  R.m.s.  vector  wind  speed  error  (m/s). 

Tables  7  and  8  list  the  final  r.m.s.  deviations  of  height 
and  wind  fields  from  observed  data,  after  five  analysis  passes  with  a 
final  radius  of  influence  R=2d.  The  reason  why  lOOOmb  height  and  wind 
errors  and  lOOmb  wind  errors  are  so  small  is  that  the  analysis  is 
constrained  by  very  few  reports  at  these  levels.  This  serves  to 
Illustrate  the  tongue-in-cheek  forecaster's  maxim,  that  nothing  can 
ruin  a  good  analysis  like  more  data. 
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Appendix  C 

Vertical  Velocity:  Formulation,  Errors,  and  Consistency 


This  appendix  describes  the  methods  for  computing  boundary 
conditions  on  omega,  the  finite-difference  forms  of  the  kinematic  and 
quasi -geostrophlc  equations,  and  such  elements  as  surface  pressure, 
divergence,  vorticity,  and  deformation.  When  calculating  horizontal 
derivatives  on  a  square  grid  in  polar  stereographic  projection,  it  is 
necessary  to  calculate  a  local  value  of  grid  distance 
l+slniji 

D=  d  x( - )  (C. 

l+s1n60° 

as  a  function  of  the  latitude,  $,  and  the  constant  grid  length 
d=190.5km  valid  at  60°N. 

The  adiabatic  method.  The  finite  difference  form  for  omega 
by  the  adiabatic  method  (assuming  dry-adiabatic  motion)  Is: 

Rd  T(t+3hr)-T(t-3hr)  g 

«j=—  {  (  -  —  ■  — . —  )  +  -  —  - (  -AyZAvT  +AXZAyT)}  (C. 

op  6  hr  f0(2D)z 

The  operators  Ax  and  av  refer  to  finite  differences  over  two  horizontal 

aT 

grid  distances.  For  the  first  and  last  data  sets,  --  over  the  nearest 

3t 

three  hours  was  assumed  to  approximate  the  value  at  the  end  times. 

The  adiabatic  method  is  sensitive  to  observational  error. 


especially  in  the  local  rate  of  change  of  temperature.  Let  the 


root-mean-square  (r.m.s.)  errors  of  the  SESAME  observations  at  lOOmb 
(Fuel berg  (1974))  be  ot-l°C  and  oz=30m.  For  fo=.92xl0"4s-1,  p=104  Pa, 
oslSOXlO^nP/kg/Pa,  2D=322km,  and  time  differencing  over  six  hours, 
one  would  expect  an  error  In  the  first  term  of  equation  (C.2)  of  .17ub/s, 
and  In  the  second  term  of  .0045ub/s.  Thus,  one  would  expect  errors  In 
lOOmb  omega  computed  by  the  adiabatic  method  as  large  as  ±.17yb/s. 

Derived  surface  quantities.  Although  surface  temperature  was 
grldded  directly  from  rawlnsonde  data,  surface  pressure  ps  and  the 
resultant  surface  density  ps  were  analyzed  indirectly  from  gridded 
values  of  the  mean  virtual  temperature  in  the  lowest  lOOmb  above  the 
surface,  Tv.  This  mean  virtual  temperature  was  assumed  to  be  more 
conservative  than  observed  surface  pressure  and  temperature,  especially 
In  the  mountains. 

PS  «  P  exp(— (Z-Zs))  (C.3) 

RdTv 

p$  =  Ps  !  ^d^s  (C.4) 

where  Z  is  the  vertically  Interpolated  height  lOOmb  above  the  surface 
Zs,  and  Ts  is  the  gridded  surface  temperature. 

Since  the  terrain  height  varied  so  much  over  the  analysis  domain, 
the  standard  approach  of  calculating  geostrophlc  surface  winds  from  a 
lOOOmb  height  field  obtained  by  downward  extrapolation  did  not  seem 
justifiable.  Instead,  a  horizontal  pressure  field,  P,  was  computed 
locally  at  each  grldpolnt  at  the  terrain  height  of  the  grldpoint. 

1  1  AyP 

Then  uq= - —  -  (C.5) 

*  2  (2D)  fops 
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where  fo=.92X10"4s-1  Is  a  constant  value  of  the  coriolis  parameter  and  P  is 
the  local  pressure  field.  To  find  P,  first  a  pressure  level  p0  was  located 
above  the  gridpoint  terrain  height,  Z0,  at  which  all  four  of  the  closest 
gridpoints  in  the  x  and  y  directions  were  above  ground.  For  each  of 
these  neighboring  points,  the  mean  virtual  temperature  in  the  layer  above 
the  surface  Is  taken  either  as  t=Tv  (discussed  above)  when  p0  is  no  more 
than  lOOmb  above  the  surface,  or  as  x=^(Tv+T),  where  T  is  the  average 
temperature  in  the  lOOmb  layer  below  p0.  Refer  to  Figure  25.  Then  for 
each  of  the  four  neighboring  points, 

P i=P0  exptJLuf-Zo)),  1*1,4  (C.7) 

RdT 

from  the  grldded  value  of  Zf  at  pressure  level  p0  at  the  respective  neigh¬ 
boring  points.  To  repeat,  Z0  is  held  constant  at  the  terrain  height  of  the 
center  gridpoint  where  tg  will  be  calculated.  These  four  values  of 
Pi.  i =1 ,4 ,  are  used  for  one  value  of  ug  and  vg  In  equations  (C.5)  and  (C.6). 

The  kinematic  method.  The  finite-difference  equation  used  to 
calculate  divergence  from  observed  wind  components  is: 

1 

6  *  1  (AvU^&yV )  (C.B) 

(2D)  y 

Let  us  assume  that  the  average  vertical  column  in  the  kinematic  calcu¬ 
lations  ran  from  800mb  to  lOOmb,  and  try  to  estimate  the  rms  error  which 
will  develop  In  the  lOOmb  omega  field  which  will  be  calculated.  Now,  the 

wind  errors  In  two  vertically  adjacent  gridpoints  should  be  positively 
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correlated,  assuming  the  errors  are  generated  consistently  during  each 
rawlnsonde  ascent.  The  positive  correlation  need  not  be  large,  only 
greater  than  zero.  Now  the  root-mean-square  error  of  any  quantity  (a+b) 
averaged  or  integrated  from  two  fields  a  and  b  is  related  to  the 
Individual  rms  errors  by: 

°a+b  3  j|j  IC«(a+b)1  )2  =  j-  paf+fib,)2 

3  l  }((5ai)2  +  (2«ai«bi)  +  (6b,)2)  ( C - 9 ) 

3  °a2  +  °b2  *  2J  {(5ai6bi) 

Where  6aj  and  5b,  are  the  Individual  errors  over  N  gridpoints,  so  that 
6(a+b)j  Is  the  individual  error  of  the  sum  of  the  fields.  Since  we 
assume  a  and  b  are  positively  correllated,  £(6aj5bi)>0,  implying  that 

°a+b>0a2+0b2  (c*10) 

For  a  column  of  seven  layers,  between  SOOmb  and  lOOmb,  the 

variance  of  the  mean  divergence  of  a  vertical  column  Is  larger  than  the 

sum  of  the  variances  contributed  by  each  layer  of  the  atmosphere. 

>j1i  <  >  <c-ll> 

where  oA2  Is  the  variance  of  the  mean  divergence  In  the  columnn,  6k  Is 

Is  the  divergence  at  pressure  levels  po=800mb,  pi=700mb,  etc.  The  factor 
1 

of  -  accomplishes  layer-averaging.  Combining  equations  (C.8),  (C.ll),  and 
(II. 9)  and  assuming  also  that  &xu  and  AyV  each  has  an  r.m.s.  error  of  2 ov. 


Recall  that  ov,  the  rms  wind  error.  Is  a  function  of  pressure.  Using 
the  values  of  ov  in  Table  8  for  1200  GMT/10  April  1979  and  1500  GMT/ 

10  April  1979  from  800  mb  through  lOOmb,  the  expected  values  of  lOOmb 
omega  errors  would  be  around  10  ub/s. 

In  contrast,  the  values  of  lOOmb  omega  computed  by  the  kinematic 
method  had  an  r.m.s.  around  7.5  ub/s  for  these  two  data  sets,  mo?.t  of 


which  is  error  (that  Is,  adiabatic  calculations  indicate  that  lOOmb  omega 

should  be  no  greater  than  lub/s  in  absolute  value).  The  source  of  this 

discrepancy  (ou=85ub/s  vs  au-7.5ub/s)  lies  in  the  assumption  that  the 

horizontal  gradients  of  u  and  v  have  r.m.s.  errors  comparable  to  2av. 

One  must  conclude  that  Axu  and  AyV  each  has  an  r.m.s.  error  no 
2 

more  than  --ov.  It  is  certainly  reasonable  to  infer  that  for  a  Cressman 
analysis,  the  horizontal  gradient  of  any  scalar  has  an  rms  error  on  the 
order  of  1/10  the  error  of  the  scalar  itself. 


O'Brien's  technique. 


factor,  k(k+1j 

*  ( k) = -■  — 
K(K+1) 


O'Brien's  (1969)  quadratic  correction 
k»0,K  (C.13) 


depends  somewhat  on  the  total  number  of  vertical  grid  spaces,  K,  but  it 
decreases  rapidly  for  all  K  to  values  less  than  34%  at  500mb  and  below 
(see  Figure  26). 
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The  quasi -geostrophic  method.  For  Indices  i»  j,  and  k  in  the 
x,  y,  and  (-p)  directions,  the  finite  difference  form  for  the 
residual  of  the  QG  omega  equation  (11.12)  is: 

r1,j,k  =  (“1+l,j,k  +  “1,j+l,k  +  w1-l,j,k  +  “1J-l,k  "  ^i.j.k) 
f  202 

+  - — ?  (“U,k+1  +  “i,j,k-l  ’  2“i  j  k) 

a(dp)2 

(C-14) 

Here  Ff  is  a  pre-defined  forcing  function,  dp  Is  always  200mb,  and  all 
values  of  omega  are  initialized  at  zero.  Using  an  overrelaxation  factor 

1  3 

a  *  2/{l+s1n  arccos  (-  l  cos  (*/N<))}  (C.15) 

3  1=1 

where  Nj  is  the  number  of  grid  intervals  in  each  direction,  one  can  make  a 
better  guess  of  omega: 

f02D2 

“  i.j.k  a  “1  ,j,k  +  n.j.k  a/(4+2“— 3)  (C.16) 

The  forcing  function  is  fully  consistent  within  itself  only  if 

(1)  second  derivatives  are  calculated  over  four  gridlengths  in  combination 

l+sin<|> 

with  first  derivatives  over  two  gridlengths,  and  (2)  the  factor  ( —  -) 

l+s1n60# 

Is  applied  to  the  reference  gridlength  d=190.5km  locally,  at  the  center 
of  each  horizontal  finite  difference  calculated.  Since  it  was  deemed 
necessary  to  violate  (1)  In  order  the  maximize  the  interior  domain  where 
values  of  Fj  j  could  be  calculated,  an  alternative  to  (2)  was  required 
as  well. 

Specifically,  laplaclans  were  calculated  over  two  grid  lengths 
,2zi.j  *  +  zt-i,j ♦  zi,j+i  *  z),m  • 
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s 


I 


I 

-» 

I 

E 


and  the  factor( - )*  was  applied  to  each  value  of  the  forcing 

l+sin60“ 

function  after  all  the  Intermediate  calculations  were  made  using  a  constant 
d=190.5  km.  In  addition,  the  second  vertical  derivative  was  computed 
Inconsistent  with  the  first  derivative: 

82u  1 


*  - 9(<*>k+l+“k-l-2"k) 

(dp)2 


(C.18) 


Since  relative  vorticity  has  not  been  corrected  for  latitude 

c  =  — v2x  (C.19) 

*o 

the  coriolls  parameter  must  have  this  correction  factor  removed: 

l+s1n$ 

f  =  14.292  sin*/( - )2  (C.20) 

l+s1n60° 

The  finite  difference  forms  of  the  four  quasi -geostrophlc 

l+sin$ 

forcing  functions  are  (dp=200mb,  Ks(  ■  ■)  ,  d=190.5km): 

l+sin60° 


Fl=—  — (af-ah)  where  (C.21) 

o  (dp) 

9  1 

a  = - -(+AyZAxU+f)-AxZAv(c  +  f)) 

f0  (2d)2  *  * 


from  the  vorticity  advectlon  above  (at)  and  below  (a^)  the  grldpolnt. 


bij=^7[7^(AyZ(Zi+l,j,k+1'Zl+1,j,lc_1‘Zl‘1’j,,c+1+Zi‘1,j’k‘1> 

-  Axz^zi,j+l,k+l"zi,j+l,k-l"zi ,j-l,k+l+21 ,j-l,k-l)} 

where  b  Is  horizontal  temperature  advection. 

9  K4 

F3  =~T7T7ZTiAxn  ( Zi .  j+1  *k+1"Zi  *  j+1  .k-l‘z1 ,  j-1  ,k+l+2i ,  j-1  ,k-l  > 

a  ( dp )  ( 2d) 41 

-Ayn(zi+1j>k+i-z1+1jtk-rzi-l,j,k+l+zi-l,j,k-l))  (c-23> 


where  nifj  *  2^  j  +  f1fj 
-g  K4  1 

F4 '  r-  tzt*  r7T(E<.j.|‘(Di.j.k-i  -  “nj.w1 

f0o  (2dT  (dp) 

’  Di,j,k  (E1,j,k-l  "  E1J,k+l))  (c-24) 

where  0  and  E  are  the  components  of  horizontal  deformation.  Since  F4  can 


be  calculated  either  as  a  residual  of  the  other  forcing  functions 

E4i,j,k  =  F11,j,k  +  E2t,j,k  ”  E31,j,k  (C.25) 

or  "independently"  as  a  separate  function  by  equation  (C.24)  it  was  not 

necessary  to  make  F4  consistent  with  the  other  forcing  functions  so  long 
as  the  resultant  omega  fields  agreed  to  a  reasonable  extent.  It  seemed 
appropriate  to  compute  shearing  deformation  In  a  way  which  is  not  con¬ 
sistent  with  derivatives  over  two  grldlengths.  The  stretching  deformation 
Is  computed  from  four  points  at  a  distance  (/2d)  from  the  center  point: 
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(C.26) 


2g 

E1,j=777T2(Zi+1,j+1_Zi+1.j-1‘Zi-1J+1+Zi-1,M) 

V2d > 

A  consistent  finite  difference  would  be  computed  from  four  points  at  a 
distance  (2d)  from  the  center: 

9 

Di,j=^^2(Zi+2J-ziij+2-21,j-2+zi-2,j)  {C*27) 

One  senses  a  scale  inconsistency  between  the  D  and  E  components,  as 
illustrated  in  Figure  27.  Considering  the  perfect  similarity  of  the 
two  operators  (C.26)  and  (C.27),  the  analyzed  fields  will  not  have  the 
same  degree  of  accuracy  unless  D  is  computed  from  four  points  at  a  dis¬ 
tance  (/2d)  from  the  center.  The  finite  difference  form  for  D  used 
in  this  research  was  based  on  linear  Interpolation 
2g 

Di*j=r^7r2(C/7”1](Zi+2'j'Zl^+2'zi*j'2^i-2,j) 

f0(2dp 

+  [2-/U(zt+1  J-*1J+l-z1,j-l+z1-l,j))  <C.28) 

Notice  that  the  factor  2  appears  in  the  numerator  analogous  to  equation 
(C.26). 
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Appendix  D 

Sanders'  Analytic  Model 


The  model .  An  Interesting  experiment  was  performed  in  which  a 
known  height  field  with  a  known  quasi -geostrophic  (QG)  omega  field  was 
compared  with  the  product  of  Cressman  analysis  and  numerical  solution  of 
the  QG  omega  equation. 

Sanders  (1971)  created  an  analytic  model  of  an  atmospheric 
disturbance  and  solved  the  QG  omega  equation  analytically  for  a  three- 
dimensional  omega  field.  A  hypothetical  vertical  temperature  profile, 
Tm(p),  is  varied  sinusoidally  in  x(east)  and  y(north)  by  the  function 

1000  2k  2n 

T(x,y,p)=Tm(p)-{l-aln - )x(Ay+Tacos— x  cos— y)  (D.l) 

P  L  L 

where  the  amplitude  of  the  oscillations  are  specified  (Ta),  and  the 
oscillations  are  damped  with  height  to  zero  at  the  tropopause.  The 
constant  a  determines  the  tropopause  pressure,  above  which  the  reversed 
temperature  perturbations  become  unrealistically  large.  The  model 
requires  a  synoptic  scale  wavelength,  L,  to  scale  both  the  north-south 
and  east-west  perturbations,  and  a  zonal  temperature  gradient,  A. 

Hydrostatic  integration  of  the  mean  temperature  profile, Tm( p ) , 
yields  a  mean  height  profile,  Zm(p).  Now  we  construct  a  surface  (lOOOmb) 
height  field  (Zm(1000)*0) : 


6 


(D.2) 


2*  2i 

Z{  x, y,  1000 ) =zacos— (x+X )  cos— y 
L  L 


where  Za  Is  the  amplitude  of  the  surface  wave  and  X  is  the  phase  difference 
between  the  surface  temperature  and  geopotential  perturbations.  The 
vertical  integral  of  temperature  gives 

2k  2k 

Z(x,y,p)=zm(p)  +  zacos— (x+x)cos— y  (D.  3) 

L  L 

1  1000  1000  2*  2w 

-  -CR^ln - (Rja^Hln - r  ](Ay+Tacos— x  cos— y) 

g  p  p  L  L 


Assuming  zero  vertical  motion  at  p=1000mb  and  p=0,  Sanders  (1971)  solved 
for  omega 

2n  2w 

w(x.y.p)=  <*>i(p)  sin— x  cos— y  (D.4) 

L  L 

2k  2k 

+-ui2 ( P )  sin— (x+X)  cos— y 
L  L 

4ir 

+u>3( p)  sin— y 
L 

where  the  three  pressure-dependent  omega  functions  have  quite  compli¬ 
cated  algebraic  expressions.  Figures  28a  and  b  illustrate  the  analytic 
height,  temperature,  omega,  and  (derived)  relative  vorticity  fields 
generated  by  equations  (0.1)  through  (D.4). 

The  analysis.  Sanders'  equations  were  used  to  create  gridded 
values  of  height  and  omega,  as  well  as  synthetic  height  soundings  at 
the  location  of  all  47  stations  of  the  SESAME  data  set  at  1200  GMT  on 
10  April  1979.  The  constants  used  in  equations  (D.1-D.4)  were  chosen 
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identical  to  "Case  C"  in  Sanders'  (1971)  article,  except  that  the  mean 
temperature  profile,  Tm{p),  was  taken  from  the  NACA  Standard  atmosphere 
(Haltlner  and  Martin  (1957)). 

Tm(p)  (NACA  Standard  temperature  profile) 
a  .722  (tropopause  at  250  mb) 

A  1.12°C  per  100  km 

Ta  7.25° 

Za  104  gpm  ~  1020m^/s^ 

L  2900  km 

X  L/4  =  725  km  (1000  mb  cold  center  west  of  1000  mb  low) 

Table  9:  Constants  supplied  to  Sanders'  analytic  model. 

The  objective  analysis  of  the  47  synthetic  soundings  based  on 
heights  alone  was  quite  poor.  The  worst  analysis  was  at  lOOOmb: 

29m  r.m.s.  deviation  of  the  gridded  height  analysis  from  the  synthetic 
lOOOmb  height  data.  But  the  height  field  analyzed  with  the  help  of 
synthetic,  quasi-geostrophic  winds  was  superb:  2m  r.m.s.  height  deviation 
at  lOOOmb.  The  wind  analysis  (from  soundings  of  geostrophic  wind  computed 
from  Sanders'  analytic  height  field  differenced  over  100km  interals)  had 
r.m.s.  errors  between  2m/s  at  the  surface  and  3.6m/s  at  300mb.  The  wind- 
assisted  height  analysis  at  500mb  is  shown  in  Figure  29a. 

Next,  two  forms  of  Sanders'  height  field  were  tested  in  the 
quasi-geostrophic  ( QG )  programs.  The  height  field  objectively  analyzed 
from  synthetic  soundings  was  used  to  compute  QG  omega  for  comparison 
against  the  analytic  QG  omega  of  equation  (0.4).  The  analytic  height 
field  of  equation  (0.3)  was  used  to  compute  the  four  separately  forced 


Q6  components  FI,  F2,  F3=(2A+C),  and  F4=(-2a).  In  all  calculations, 
top  and  bottom  boundary  conditions  were  taken  as  zero,  exactly  as  in 
Sanders'  model. 

Figure  29b  shows  how  well  the  500mb  QG  omega  calculation 
approximates  the  analytic  QG  field,  except  close  to  the  side  boundaries, 
where  omega  was  held  to  zero.  The  partitions  of  QG  omega  in  Figure  30a-d 
are  all  reasonable.  The  most  notable  feature  is  the  magnitude  of  vertical 
motion  due  to  F4,  Winn-Niel sen' s  deformation  function,  which  is  at  least 
an  order  of  magnitude  less  than  the  other  three  components  of  QG  omega. 
The  location  of  this  F4- induced  omega  relative  to  the  synoptic  situation 
is  also  of  interest. 
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Fiq.lt  Computational  grid,  SESAME  ravin-  Fig. 2 5  Time -averaged  700  mb  height 

•onde  stations  (not  underlined),  and  supple-  10-11  April  1979.  Subjective  corrections 

mental  ravinsonde  stations  (underlined).  shown  as  dashed  lines. 


Fig. 5a:  Surface  fronts  and  Fig. 5b:  Nephanalysis  at  1200 

alt.  setting  (102in  Hg:  solid)  GMT  10  April  1979  from  GOES  IR 

and  quasi -geos trophic  surface  data,  with  convective  and  sandstorm 

omega  (yb/s:  dashed)  at  1200  GMT  activity  from  radar  and  hourlies. 
io  April  1979.  Ceilings  (scalloped). 


u 


V  / 


>£< 


,  500  ^  height  („: 

dLi^llneai'  te'^ratur«  <C:  heavy 
/,n-Sd-1  and  relatlv®  vorticity 

*  s  light  dashed)  at  1200  GMT 


Minor  wave  trough  ( 


\  71  < 


Fig.  Jds  500  mb  quasi- 

geostrophic  omega  (pb/s)  at 
1200  GMT  10  April  1979. 


c 


Fig,5e»  500mb  kinematically 
computed  omega  U  b/s)  corrected  by 
O’Brien’s  method,  at  1200  GMT 
10  April  1979. 
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Fig.  :  SOOmb  omega  (ub/s) 

due  to  the  laplacian  of 
temperature  advection  (F2),  at 
1200  GMT  10  April  1979. 


Fig.5h:  SOOmb  omega  (wb/s) 

due  to  the  advection  of  vorticity 
(2C^f)  by  the  thermal  wind  (F3), 
at  1200  GMT  10  April  1979. 


Fig.5ij  SOOmb  omega  (wb/s) 

due  to  Wiin-Nielsen's  deformation 
function  (F4)f  at  1200  GMT 
10  April  1979. 


Fig.5f«  500mb  omega  (pb/s) 

due  to  differential  vorticity 
advection  (FI),  at  1200  GMT 
10  April  1979. 


Pig.  6a  «  Surface  fronts  and  Pig. 6b:  Nephanalysis  at  1500 

alt.  setting  (102in  Hg:  solid)  GMT  10  April  1979  from  GOES  video 
and  quasi -geos trophic  surface  data,  with  convective  and  sandstorm 
omega  (ub/s:  dashed)  at  1500  Q4T  activity  from  radar  and  hourlies. 
io  April  1979.  Dry  line  \,Cei lings  (scalloped). 


rwi 


rig.OCt  500  mb  height  (m: 

solid  lines),  temperature  (Cs  heavy 
dashed),  and  relative  vorticity 
<10  s  s  liqht  dashed)  at  ihoo  GMT 
Minor  wave  trough  (  .  ) . 


Fig.6di  500  mb  quasi- 
geostrophic  omega  (ub/s)  at 
1500  GMT  10  April  1979. 


rig. 500nb  kinematically 
computed  omega  (wb/s)  corrected  by 
O'Brien’s  method,  at  1500  GMT 
10  April  1979. 


Fig.  7a  i  Surface  fronts  and  Fig.?b»  Nephanalysis  at  1800 

alt.  setting  (I02in  Hg:  solid)  GMT  10  April  1979  from  GOES  video 
and  qcasi-geostrophic  surface  data,  with  convective  and  sandstorm 
omega  (wb/s:  dashed)  at  1800  GMT  activity  from  radar  and  hourlies. 

10  April  1979.  Ory  line  XCeilings  (scalloped). 
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Pig. 8a  |  Surface  fronts  and 
alt.  setting  (lO^in  Hg:  solid) 
and  quasi -geostrophic  surface 
omega  (yb/s:  dashed)  at  2100  GMT 
10  April  1979.  Dry  line  ( 


Fig. 8b i  Nephanalysis  at  2100 

GMT  10  April  1979  from  GOES  video 
data,  with  co  "'active  and  sandstorm 
activity  from  radar  and  hourlies. 

).Ceilings  (scalloped)  and  cold¬ 
est  clouds  (solid  outline). 


Fig.Qci  500  mb  height  (m: 

■olid  lines),  temperature  (C:  heavy 
dashed),  and  relative  vorticity 
(10  s  liqht  dashed)  at  2100  GMT, 
Minor  wave  trough  (  ) . 


fig,  Q0  ,  SOOmb  kinematically 

computed  omega  (wb/s)  corrected  by 
O’Brien's  method,  at  2100  GMT 
10  April  1979. 
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Fig  «8d  j  500  mb  quasi- 

geostrophic  omega  (i>b/s)  at 
2100  GMT  10  April  1979. 
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Fig  .©f »  500mb  omega  (ub/s) 

due  to  differential  vorticity 
advection  (FI),  at  2100  GMT 
10  April  1979. 


Fig.Bg:  500mb  omega  (wb/s) 

due  to  the  laplacian  of 
temperature  advection  (F2),  at 
2100  GMT  10  April  1979. 
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500mb  omega  (i  b/») 
due  to  the  advection  of  vortici  ty 
(2C«-f)  by  the  thermal  wind  (F3), 
at  2100  GMT  10  April  1979. 


Fig. Si  i  oOOmb  omega  Cub, •) 

due  to  Wi in-Nie 1  sen ' s  def srma * i on 
function  (F4),  at  2^00  SMT 
10  April  1979. 


Pig. 9a  I  Surface  fronts  and  Fig.  9b  t  Nephanalysis  at  0000 

alt.  setting  (102in  Hg:  solid)  GMT  11  April  1979  from  GOES  video 

and  quasi -geostrophic  surface  data,  with  convective  and  sandstorm 
omega  (pb/ss  dashed)  at  0000  GMT  activity  from  radar  and  hour lies. 

11  April  1979.  Dry  line  ( ).  Ceilings  (scalloped)  and  cold¬ 
est  clouds  (solid  outline). 


Pig. 9c  t  50 0  mb  height  (ms 

•olid  lines),  temperature  (C:  heavy 
dashed),  and  relative  vorticity 
<10“5s”1s  light  dashed)  at  0000  GMT. 
Minor  wave  trough  (  )• 


rig.9d»  500  mb  quasi” 
geostrophic  omega  (Wb/s)  at 
0000  GMT  11  April  1979. 


Fiq.9et  500ab  kinematically 

computed  omega  (ub/a)  corrected  by 
O'Brien’ a  method,  at  0000  GMT 
11  April  1979. 


rig. 1  Oat  Surface  fronts  and  rig. 1  Obi  Nephanalysis  at  0300 
alt.  setting  (102in  Hg:  solid)  «T  11  April  1979  from  GOES  IR 
and  quasi -geos trophic  surface  data,  with  convective  and  sandstorm 
omega  (wb/s:  dashed)  at  0300  GMT  activity  from  radar  and  hourlies. 

11  April  1979.  Dry  line  l^^).Ceilings  (scalloped)  and  cold¬ 
est  clouds  (solid  outline). 
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Fig.lOd*  WO  »b  qu*8i* 
geo*trophie  omega  (ub/s)  at 
C500  GMT  11  April  1979. 


Fig.lOei  SOOmb  kinematically 
computed  omega  (yb/a)  corrected  by 
O'Brien’*  method,  at  0900  GMT 
11  April  1979. 


c 


Fig.lla*  Surface  fronts  and  Fig.il bt  Nephanalysis  at  0600 

•It.  setting  (102in  Hg:  solid)  GMT  11  April  1979  from  GOES  IR 

and  quasi -geos trophic  surface  data,  with  convective  and  sandstorm 
omega  (ub/s:  dashed)  at  0600  GMT  activity  from  radar  and  hourlies. 

11  April  1979.  Dry  line  ).  Ceilings  (scalloped)  and  cold¬ 

est  clouds  (solid  outline). 


Fig .11c  i  500  ab  height  (m: 

solid  lines)*  temperature  CCs  heavy 
dashed),  and  relative  vorticity 
(10"5s_1:  light  dashed)  at  0600  GMT. 
Minor  wave  trough  <  .A  4l—  )  • 


Fig  .lid  i  500  mb  quasi- 

geostrophic  omega  (Wb/s)  at 
0600  GMT  11  April  1979. 


Fig.llBi  500ab  kinematically 
computed  omega  (pb/s)  corrected  by 
O'Brien's  method,  at  0600  GMT 
11  April  1979. 


Pig. 12a:  Surface  fronts  and 

•It.  setting  {102in  Hg:  solid) 
•nd  quasi -geos trophic  surface 
omega  (pb/s:  dashed)  at  0900  GMT 
**  April  1979.  Dry  line  ( 


Pig. 12b:  Nephanalysis  at  0900 

GMT  11  April  1979  from  GOES  IR 
date,  with  convective  and  sandstorm 
activity  from  radar  and  hour lies. 

). Coldest  clouds  (solid  line). 


Fig.  12 Cl  500  ab  height  (a: 

^Kd^inea)'  tem?«rature  (C:  heavy 
das  and  relative  vortieity 
(10  »"  s  light  dashed)  at  0900  GHT, 
Minor  wave  trough  (  . 1 . 


Pig.l2di  500  mb  quasi-  rig.l2e«  500mb  kinematically 
geostrophic  omega  (wb/s)  at  computed  omega  (ub/s)  corrected  by 
0900  GMT  11  April  1979.  O'Brien’s  method,  at  0900  GMT 

11  April  1979. 
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Pig. 13a:  Surface  fronts  and  Pig*13b:  Nephanalysis  at  1200 

alt.  setting  (102in  Hg:  solid)  GHT  11  April  1979  from  GOES  IR 
and  quasi -geos trophic  surface  data,  with  convective  and  sandstorm 
omega  (iib/s:  dashed)  at  1200  GMT  activity  from  radar  and  hourlies. 

11  April  1979.  Dry  line  KColdest  clouds  (solid  line). 


Fig.l3f»  500mb  omega  (Wb/e)  Fig.l3g*  SOOmb  omega  (ub/s) 

due  to  differential  vorticity  due  to  the  laplacian  of 

advection  (FI),  at  1200  GMT  temperature  advection  (F2),  at 

11  April  1979.  1200  GMT  11  April  1979. 


Fig.  13^1  SOOmb  omega  (ib/i)  Fig.l3ii  SOOmb  omega  (wb/e) 
due  to  the  advection  of  vorticity  due  to  Viin-Nielsen's  deformation 
(2C+f )  by  the  thermal  wind  (F3),  function  (F4),  At  1200  GMT 
at  1200  GMT  11  April  1979.  11  April  1979. 
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Fig.l4ai  Richardson  numbers 
(see  text).  Regions  at  and 
above  500mb  at  1 200  GMT  1 0 
April  1979  with  0*Ri<.95 
(solid)  and  0<Rj< .25 (dashed ) . 
500mb  trough  (  )  . 


Fig.l4b«  Regions  below  500mb 
at  1200  GMT  10  April  1979 
with  (XRi*.95(solid)  and 
0-iRj^.25tdashed) .  Surface 
Ri<*95(scalloped) .  Surface 
fronts  and  lows. 


Fig.l5ai  Richardson  numbers 
(see  text).  Regions  at  and 
above  500mb  at  1500  GMT  10 
April  1979  with  0<Ri^.95 
(solid)  and  O^^Sfdashed) 
500mb  trough  (  )  . 


at  1500  GMT  10  April  1979 
with  0*Rj>.95(solid)  and 
0<Ri< .25( dashed) .  Surface 
Ri^.95(scallaped) .  Surface 
fronts,  low,  and  dry  line. 
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Fig. 16a*  Richardson  numbers 
(see  text).  Regions  at  and 
above  500mb  at  1  S00  GMT  10 
April  1979  with  0*Ri<.95 
(solid)  and  0<R^ .2$( dashed ) . 
500mb  trough  (  ) , 


a  Jr 

Fig»1?a«  Richardson  numbers 
(see  text).  Regions  at  and 
above  500mb  at  21 00  GMT  10 
April  i 979  with  0*RjS.95 
(solid)  and  0^^^ .25 (dashed ) . 
500mb  trough  (  ) . 


Fig. 1 6b 1  Regions  below  500mb 
at  1300  GMT  10  April  1979 
with  0<£R<  *.95(solid)  and 
0*Rj^.25tdashed) .  Surface 
.95(scalloped) .  Surface 
fronts,  low,  and  dry  line. 


Fig.l7b»  Regions  below  500mb 
at  2100  GMT  10  April T979 
with  0*R<  $.95(solid)  and 

25* dashed ) .  Surface 
Ri^»95(scalloped) .  Surface 
fronts,  lows,  and  dry  line. 


Fig. 1 8a*  Richardson  numbers  Fig. 19a*  Richardson  numbers 

(see  text).  Regions  at  and  (see  text).  Regions  at  and 

above  500mb  at  0000  GMT  11  above  500mb  at  0300  GMT " 1 1 

April  1979  with  O^R?<.95  April  1979  with  0<Ri<.95 

(solid)  and  0<Ri£.25(dashed) .  (solid)  and  0<Rj£.25(dashed) 
500mb  trough  (  ) .  500mb  trough  ( ) . 


Fig.l8bi  Regions  below  500mb  Fig. 19b*  Regions  below  500mb 
at  0000  GMT  11  April  1979  at  0300  GMT  11  April  1979 

with  0<RjS.95(solid)  and  with  0<R<* .95(solid)  and 

0<Rj*.25(dashed) .  Surface  0<Ri$.25( dashed) •  Surface 

Ri^»95(scalloped) .  Surface  Rj£.95(scalloped) .  Surface 

fronts,  lows,  and  dry  line.  fronts,  lows,  and  dry  line. 
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Fig. 20a :  Richardson  numbers 
(see  text).  Regions  at  and 
above  500rab  at  0600  GMT  11 
April  1979  with  0  R<  .95 
(solid)  and  Q^R^. 2$ (dashed ) . 
500mb  trough  ( ) , 


Fig. 21a i  Richardson  numbers 
(see  text).  Regions  at  and 
above  500mb  at  0900  GMT  11 
April  1979  with  0<R;£.95 
(solid)  and  0<Rj£. 2$ ( dashed ) . 
500mb  trough  (  ) . 


gggj 

H 


~  Y\  T-  s' 


Fig. 22a i  Richardson  numbers 
(see  text).  Regions  at  and 
above  500mb  at  1200  GMT  11 
April  1979  with  O^R^.95 
(solid)  and  0*Rj<. 25 ( dashed ) . 
500mb  trough  ( ) . 


Fig. 22b i  Regions  below  500mb 
at  1200  GMT  11  April  1979 
with  0<R;£.95(solid)  and 
0<Ri< *25( dashed) .  Surface 
•95(scalloped) .  Surface 
fronts,  lows*  and  dry  line. 
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Pig. 23 >  Cressman's  weighting  function  versus  the  ratio 
of  distance  froa  the  gridpoint  (d)  to  the  radius  of 
influence  (R). 


c4*  *3* 

** . i--. 

.I  * 


Rig.24i  Bilinear  interpolation  of  the  value  of  a  scalar 
to  a  point,  p0,  froa  gridded  values  at  the  four  closest 
gridpoints,  gj  through 
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FigOOai  SOOmb  omega  (ub/s>  Fig.30bi  SOOmb  omega  (ub/s) 

due  to  differential  vorticity  due  to  the  laplacian  of  thermal 

advection  (Fi)  for  Sanders'  analytic  advection  (F2)  for  Sanders'  analytic 
height  field  (Figure  E.la).  height  field  (Figure  E.la). 
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Fig. 30c  «  SOOmb  omega  (gb/s) 

A**  to  the  advection  of  vorticity 
J  by  the  thermal  wind  (F3) 
for  landers'  analytic  heights. 


function  (F4)  for  Sanders'  analytic 
height  field  (Figure  E.la). 
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